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The large-distance dynamics in quarkonium systems is investigated, in the large N c limit, 

through the saturation of Wilson loop averages by minimal surfaces. Using a represen- 
ts ■ 

q ■ tation for the quark propagator in the presence of the external gluon field based on the 

use of path-ordered phase factors, a covariant three-dimensional bound state equation of 
the Breit-Salpeter type is derived, in which the interaction potentials are provided by 

Q-i' 

qj \ the energy-momentum vector of the straight segment joining the quark to the antiquark 

and carrying a constant linear energy density, equal to the string tension. The interac- 
tion potentials are confining and reduce to the linear vector potential in the static case 
and receive, for moving quarks, contributions from the moments of inertia of the straight 
segment. The self-energy parts of the quark propagators induce spontaneous breakdown 
of chiral symmetry with a mechanism identical to that of the exchange of one Coulomb- 
gluon. In the nonrelativistic limit, long range spin-spin potentials are absent; the moments 
of inertia of the straight segment provide negative contributions to the spin-orbit poten- 
tials going in the opposite direction to those of the pure timelike vector potential. In the 
ultrarelativistic limit, the mass spectrum displays linear Regge trajectories with slopes in 
agreement with their classical relationship with the string tension. 
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1 Introduction 



The Wilson loop £Q appears as one of the most efficient tools for probing the large- 
distance properties of QCD. It provides a natural criterion for confinement through the 
area law and also participates as a basic ingredient in the formulation of lattice gauge 
theories PjJ |2] . Equations concerning path-ordered phase factors were first obtained by 
Mandelstam [2j and analyzed in the QCD case by Nambu j3]. Loop equations were 
obtained by Polyakov [Sj and Makeenko and Migdal [Bl Ej • (For reviews, see Refs. jH] and 
jn].) The loop equations actually represent an infinite chain of coupled equations relating 
vacuum expectation values of loop operators having as supports different numbers of 
closed contours emerging from contours with self-intersections. These equations must in 
addition be supplemented with constraint equations [TQJ ITT] , the so-called Mandelstam 
constraints, which are sensitive to the gauge group structure of the theory. Due to their 
complexity, the loop equations have not yet allowed a systematic resolution of QCD in 
terms of loop variables. In the two-dimensional case however, explicit expressions of the 
Wilson loop averages for various types of contour have been obtained for U(N C ) gauge 
theories ^2] • Renormalization properties of the Wilson loop averages were studied in the 
framework of perturbation theory in Refs. ^3 E] , where it was shown that the latter 
are multiplicatively renormalizable. 

Considerable simplification is obtained in the large- N c limit ^3] , corresponding to the 
planar diagram approximation of the theory. In that case, apart from the disappearance 
of many nonleading terms, it is the factorization property of the Wilson loop average 
for two disjoint contours that becomes mostly relevant. Makeenko and Migdal studied 
the resulting equations in the above limit pointing out their equivalence with a chain of 
Schwinger-Dyson type equations O Ej . They showed that for large contours asymptotic 
solutions exist correponding to the minimal surfaces bounded by the loops. 

Concerning physical applications, Eichten and Feinberg, using the area law for the Wil- 
son loop, could obtain the general expression of spin- dependent forces for quark-antiquark 
systems to order 1/c 2 in terms of color electric and magnetic field correlators ^Hj. This 
problem was also investigated by Gromes [17j . Later on, Prosperi, Brambilla et al. com- 
pleted these results by also obtaining the velocity- dependent forces jTHl EH I2H1 E] • On 
the other hand, a wide program of investigations was undertaken by Dosch, Simonov et 
al. in the framework of the "stochastic vacuum model" [22] • 

In recent years, the Wilson loop gained renewal of interest in connection with dual- 
ity properties of different field theories in different dimensions manifested through the 
AdS /CFT correspondence [2H]. 

The purpose of the present paper is to investigate the properties of the Wilson loop 
concerning the bound state problem in QCD in the large-distance regime when the latter 



2 



is probed by minimal surfaces. Among possible nonperturbative solutions of the loop 
equations at large distances, minimal surfaces appear as the most natural ones [HI H] ; they 
produce in a simple way the area law needed for confinement and satisfy the factorization 
property, valid at large N c . A complete solution of the loop equations should necessarily 
include short-distance effects for which minimal surfaces do not seem sufficient alone to 
provide the appropriate behavior. Therefore, for that regime perturbation theory results 
should appropriately be incorporated in the Wilson loop solution; this aspect of the 
problem will not, however, be analyzed in the present paper. Other nonperturbative 
properties of the theory, related to its possible string-like behavior, might also arise from 
contributions of fluctuations of surfaces around the minimal surface bounded by the loop 

ei mum. 

The main ingredient in our approach is a representation of the quark propagator in 
the presence of an external gluon field as a series of terms involving free propagators, 
path-ordered phase factors along straight lines and their derivatives. That representation 
generalizes to the relativistic case the one used in the nonrelativistic limit ^H]- The gauge 
invariant quark- ant iquark Green function is then represented by a series of terms involving 
Wilson loops having as contours skew polygons with an increasing number of sides. Each 
Wilson loop average is then replaced by the contribution of the corresponding minimal 
surface bounded by the loop. Contrary to the usual two-particle Green functions, the 
gauge invariant Green function does not manifestly satisfy a genuine integral equation 
which might, as in the Bethe-Salpeter equation case [2Z|, result in a bound state integral 
equation. In the present case, the Green function satisfies with the Dirac operators 
two independent and compatible equations. Selecting in the large time limit the total 
momentum of a bound state [21] and taking in the center-of-mass frame the equal-time 
limit for the two particles allow, with certain mathematical assumptions, the grouping of 
terms into a form that leads to a three-dimensional Breit-Salpeter type wave equation 
pnjlHOj. where the potentials are represented by the components of the energy-momentum 
vector of the straight segment joining the quark to the antiquark and carrying a constant 
linear energy density, equal to the string tension; they involve, apart from the usual 
confining linear potential, contributions coming from the moments of inertia of the above 
segment, which plays the role of the color flux tube of the quarkonium system. 

The self-energy parts of the quark propagators, extracted from the interaction terms, 
allow the analysis of the chiral symmetry properties of the system. The situation here is 
very similar to that resulting from the exchange of one Coulomb-gluon and the latter has 
been studied in the literature. The Schwinger-Dyson equation satisfied by the self-energy 
part has a non-trivial solution leading to a spontaneous breakdown of chiral symmetry. 

From the nonrelativistic limit of the wave equation one determines the hamiltonian 
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of the system to order 1/c 2 . Spin-spin potentials are absent from the hamiltonian. The 
contributions of the color flux tube result in new terms for the orbital angular momentum 
and for the spin-orbit potentials. In particular, the momentum of the flux tube contributes 
with a negative sign to the spin-orbit potential, in opposite direction to the contributions of 
the pure timelike vector potential, and may account for the phenomenological observations 
made for fine splitting. 

At high energies, the mass spectrum displays linear Regge trajectories, the slopes of 
which tend to satisfy the classical relationship with the string tension. 

The paper is organized as follows. In Sec. 2, we review the equations satisfied by the 
Wilson loop averages and outline the particular status of minimal surfaces within the set 
of possible solutions. In Sec. 3, we display some basic properties of minimal surfaces. In 
Sec. 4, the representation of the quark propagator in the presence of an external gluon 
field is constructed. Section 5 is devoted to the construction of the representation of the 
gauge invariant two-particle Green function in terms of a series of Wilson loops having 
as boundaries skew polygons. The bound state equation is derived in Sec. 6. Section 7 
is devoted to the extraction from the interaction terms of the self-energy parts needed 
for the quark propagators. Chiral symmetry breaking is studied in Sec. 8. In Sec. 9, 
the main qualitative properties of the bound state spectrum are displayed. Summary and 
concluding remarks follow in Sec. 10. In appendix A, the main properties of minimal 
surfaces are presented. Appendix B is devoted to the determination of the normalization 
condition of the wave function. Appendix C gives details about the Breit approximation 
used for the resolution of the bound state spectrum. 



2 Loop equations 

The starting point is the gauge covariant path-ordered phase factor along a line C yx joining 
point x to point 

rv 

-ig / dz^AJz) 

U(C yx ,y,x) = U(y,x) = Pe J* , (2.1) 

where = J2 a A^t 11 , A't (a = 1, ... , — 1) being the gluon fields and t a the generators 
of the gauge group SU(N C ) in the fundamental representation, with the normalization 
trt a t b = \5 ab ■ A more detailed definition of U is given by the series expansion in the cou- 
pling constant g; all equations involving U can be obtained from the latter. Parametriz- 
ing the line C with a parameter a, C — {x(a)}, < a < 1, such that x(0) = x and 
x(l) — y, a variation of C induces the following variation of U [U(x(a),x(a')) = U(a,a'), 
^Formulas of this section are written in Minkowski space. 
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A(x(a)) = A(a)]: 



5U (1,0) = -ig5x a (l)A a (l)U(l,0) + igU(l,0)A a (0)6x a (0) 

+ig f daU(l,a)x' /3 (a)Fp a (a)5x a (a)U(a,0), (2.2) 



o 



where x' = || and F is the field strength, F^ v = d^A„ — dyA^ + iglA^ A v \. The functional 
derivative of U with respect to x(a) (0 < o < 1) is then [5]: 

= fcZ7(l, ^/(^^(^[/((T, 0). (2.3) 

Defining the ordinary derivations with the prescriptions [3] 

d r a + £ 5 

lim / da' (2.4) 



dx^(a) e^oJ a - E 5x^(a') 

lim/ da'(a'-a) . (2.5) 



one obtains: 



dx'P(a) 5x a (a) 
A derivation with respect to x fl (a) then yields 

d d <5I7(1, 0) 



^WU.^aWM)- (2-6) 



dx*(<j)dx!P{a) 5x a (a) 



igU(l,(r)(V lt Fp tx ((T))U((T,0), (2.7) 



where V is the covariant derivative, (VF) = (dF) + ig[A, F]. 

The Wilson loop, denoted $(C), is defined as the trace in color space of the path- 
ordered phase factor (j2.1|) along a closed contour C: 

1 -iq & dx^AJx) 
$(C) = — trPe "Jo , (2.8) 

where the factor l/iV c has been put for normalization. It is a gauge invariant quantity. 
Its vacuum expectation value is denoted W(C): 

W(C) = ($(C)) A , W{C 1} C 2 ) = ($(d)$(C 2 )) A , (2.9) 

the averaging being defined in the path integral formalism. More generally, one meets 
insertions of local operators 0{x) into the Wilson loop. Their vacuum expectation values 
are: 

1 , -iq <p dz^AJz 
'— tr Pp. Jc M 



(0(x)) w =(-trP(e »Jc J 0(x)))^ 



(2.10) 
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The Wilson loop and its average also satisfy equations analogous to Eqs. ()2.3|) . (J2.6|) 
and ()2.7|) . Considering in Eq. (|2.7|) three different indices and taking their cyclic per- 
mutations one obtains in the right-hand side the Bianchi identity. For the Wilson loop 
average the equation takes the form 

£ ^ a _d__d_swM = o (2 11) 

dx^dx'P 5x a 

Contraction of the indices [i and (3 in Eq. (J2.7j) leads in the right-hand side to the 
equation of motion of the gluon field. In the large N c limit the quark current can be 
neglected and the corresponding term becomes equivalent to the functional derivative 
— S/SA a which acts now on the Wilson loop. [Gauge-fixing and ghost terms mutually 
cancel each other in gauge invariant quantities f3J E] and hence can be ignored.] One 
finds: 

8 d 5W(C) = .£?N C i j_, aS 4^. _^ JTf „ ^ 1 



dx 13 dx'p 5x 



j c dy a 5\y - x) [w(C yxj C xy ) - j^W(C)} . (2.12) 



Because of the delta function, the points y that contribute to the integral are those that 
lie in the vicinity of the point x. Two cases may emerge. First, the contour C may have 
self-intersection points and point x may be one of them. In that case, there is a point y 
which coincides with x in space but is not the same one on the contour. The contours 
C yx and C xy are the complementary contours separated by the intersection point and 
contribute as independent closed contours to the Wilson loop. Second, there is always 
the contribution of the points y which are close to x in space and on the contour. In that 
case, one of the contours C yx or C xy shrinks to zero and the other one becomes identical 
to C. In the large N c limit, the Wilson loop average of the the two contours C yx and C xy 
factorizes into two loop averages and Eq. ()2.12|) becomes: 

d d 5W{C) .g 2 N c 



i dy a 5\y - x)W(C yx )W(C xy ). (2.13) 



dx 13 dx'p 5x a 2 Jc 

[The product g 2 N c is maintained fixed and finite in the above limit [TH] .] 

Equations (|2.11|) and (|2.13|) can be considered as basic equations of QCD in the large 
N c limit in loop space. The factorized structure of the right-hand side of Eq. ()2.13l) 
puts severe restrictions on the class of its possible solutions. Potential solutions of that 
equation with contours without self-intersections may not be solutions for contours with 
self-intersections. On the other hand, many physical applications related to meson spec- 
troscopy do not require consideration of contours with self-intersections (in the static 
limit). Due to the huge complexity of the treatment of the problem with general types of 
contour, it seems reasonable, in view of physical applications, to limit, at a first stage, the 
investigations to contours not having self-intersections. We shall adopt this limitation in 
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the present paper, indicating at the end of this section the way of incorporating contours 
with self-intersections. Equation (j2.13j) then becomes: 



The points y that contribute to the integral are those that lie in the vicinity of the point 
x in space and on the contour. 

A first class of solution to Eqs. (|2.11|) and (|2.13|) is provided by perturbation theory 
[T3*l ll4|. In the right-hand side of Eq. (|2.13|) . g 2 N c represents the unrenormalized coupling 
constant and the resolution of the equation is accompanied by the renormalization of the 
theory and of the Wilson loop itself. However, perturbation theory becomes unstable at 
large distances and the search for nonperturbative solutions becomes necessary. Among 
these, minimal surfaces appear, for several reasons, as natural candidates. First, they 
easily reproduce the area law related to confinement £Q ; second, they satisfy, with certain 
restrictions, the factorization law in the large N c limit; third, they are connected to the 
classical action of string theory, which in turn is expected to have an implicit relationship 
with QCD |U |2U 123 EH] • In the following, we concentrate on the contributions of surfaces 
to Eqs. dmH) and 

Let S be a surface of a given type, that is satisfying a given equation, and having as 
contour the closed loop C. Let A be its area. Then a representation of the Wilson loop 
average can be given by the following expression: 



where a is a constant, which will be identified with the string tension. Let us consider 
Eq. (j2.11|) . It represents the action of some local deformation of the contour at point 
x. However, any local deformation on the contour also induces an internal deformation 
of the surface itself, because the latter is constrained by its defining equation and by 
the boundary condition. It turns out that the Bianchi identity operator, present in Eq. 
(J2.11J) . is represented in the internal part of the surface (and on its contour) by the 
defining operator of minimal surfaces (surfaces having minimal area). Therefore, Eq. 
(j2.11|) implies that among surfaces of a given type, only minimal surfaces can be solutions 
of it. [The details of this assertion will be presented in the Sec. El] This result considerably 
reduces the class of surface type solutions to the loop equations. In order to incorporate 
other types of surface in representation (j2.15|) there remains the possibility of considering 
contributions of an infinite sum of all possible surfaces, which might, through mutual 
cancellations, satisfy Eq. (j2.11j) . This latter possibility was considered in Ref. [Zj. Such 
a sum can also be considered as representing fluctuations of surfaces around the minimal 




(2.14) 



W(C) = e 



iaA(S, C) 



(2.15) 
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surface and corresponding to string-like contributions [2U 123 EH] • Henceforth, we restrict 
ourselves to the study of contributions to Eq. (|2.15|) coming from minimal surfaces. 

To study Eq. ()2.14jl . one must introduce a short- distance regulator a since the right- 
hand side displays a singular behavior. Using for W(C) representation (|2.15j) with a 
minimal surface one finds that Eq. (j2.14j) is satisfied provided one has the identification 

° = (2-16) 

where C is a numerical constant. Notice that the use of an arbitrary trial functional for 
W(C) would produce in general for the right-hand side of Eq. (j2.14j) a multiplicative 
functional depending on the whole contour C, rather than a local factor depending only 
on x and on its vicinity. If the unrenormalized coupling constant g 2 N c vanished as 0(a 2 ) 
with vanishing a, then minimal surfaces would define a theory by their own with a finite 
dimensional coupling constant a. However, we know from the short- distance behavior 
of QCD, given by the perturbative solution, that the unrenormalized coupling constant 
vanishes only logarithmically. Therefore condition (j2.16|) cannot be satisfied in general. 
One must then interpret the minimal surface contribution to the Wilson loop average as 
a part of a general solution in which it represents its large-distance behavior, while the 
other part includes the perturbation theory contribution representing its short- distance 
behavior. A matching condition between the two parts at intermediate distances should 
then fix the value of the string tension a in terms of the parameter A of QCD. We shall not 
study this aspect of the problem in the present paper, but rather shall assume that such 
a matching condition exists and is fulfilled and shall henceforth consider the properties of 
the large-distance behavior of the theory as deduced from the minimal surface contribution 
to the Wilson loop average. 

Finally, let us comment on the factorization property of minimal surfaces. The minimal 
surface of two disconnected closed contours lying sufficiently far from each other is the sum 
of the minimal surfaces of each contour. This ensures the factorization property of the 
Wilson loop average valid at large N c . However, when the two contours are close to each 
other new global solutions may arise that do not lead to factorization. Thus, for instance, 
if the two contours are similar to each other and lie closely one above the other, forming 
the two bases of a cylinder, then the global (or absolute) minimal surface corresponding 
to the two contours is the area of the lateral surface of the cylinder rather than the sum 
of the areas of its two bases, a result that manifestly does not lead to factorization. In 
such cases and even for single contours having complicated forms, the minimal surface 
solution should be chosen locally (therefore not necessarily the absolute one) with global 
properties being in agreement with factorization or with other physical requirements. 

Extension of that prescription to the case of contours with self-intersections allows 
minimal surfaces to also satisfy the factorization property of the right-hand side of Eq. 
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(J2.13|) . However, when such contours lie on overlapping surfaces additional contribu- 
tions might be necessary to consider, since internal contours then belong to two adjacent 
surfaces. This is in particular the case in two-dimensional QCD |12j . 

In the next section, we shall study properties of minimal surfaces, showing some of 
the results mentioned above and deriving results needed for later investigations. 

3 Loop equations and minimal surfaces 

Let C be a closed contour and S a surface bounded by it. We shall generally consider as 
prototypes of contour those having four distinct sides, with junction points designated by 
xi, x[, x' 2 , x 2 (see Fig. Q). 




Figure 1: Prototype of a closed contour with four distinct sides. 



The surface bounded by this contour will be parametrized with two parameters, a and 
r, and a point belonging to it will be represented as x(a, r) or y(a, r). We adopt for the 
partial derivatives the usual notations 

. dx dx 
X= Ya X= d~r (3 ' 1} 

In general, the four sides of the contour will be parametrized as follows: {^iO^} with 
<T = 0, {x[x' 2 } with r = 1, {x' 2 X2} with a = 1 and {22^1} with r = 0. The area of the 
surface is 1 : 

A(C,S) = £ da £ dr(x' 2 x 2 - {x' .xf) ^ . (3.2) 

To avoid the occurrence of possible divergences when dealing with functional deriva- 
tives of the area, it is necessary to introduce a short-distance regulator in the expression 



^Formulas related to surfaces will in general be written in euclidean space. 
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of the area. We adopt for the latter the following expression, already suggested in a more 
general form in Ref. jjj: 

A= \ fda^(x) f da^(y)F(x-y), 



2 

da^(x) = dadr(x^x u - afx"), 

^ a 2 1 . . 

F = -l ^ ( 3 - 3 ) 

' (Or - yf + a?) 
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where a is the (positive) regulator and goes to zero at the end of calculations. In the limit 
a — > 0, F is actually equal to a two-dimensional 5-function: 



lim F = (x l2 x 2 - (x'.x) 2 ) 1/2 5(a - a')8(r - t') . (3.4) 

This implies that it is only points x and y lying close to each other that contribute to the 
integral and therefore a limited number of terms of the expansion of y about x have to 
be considered in general. For example, one finds the following limits: 



— I da'dr' — 



(y(a',r')-x(a,r)Y + a^' 



(x' 2 x 2 — (x'.x) 2 ) ^ , for x inside 5", 



|(x /2 x 2 — (x'.x) 2 ) ^ , for x on C. 

(3.5) 

We are interested in variations of the area of the surface S when local variations are 
introduced on its contour C . S being a surface of a given type, that is satisfying a defining 
equation, any deformation of the contour introduces corresponding deformations inside 
the surface. Let us for definiteness consider deformations on the line r = 0. The general 
deformation of the area is: 

SA = J da6x a (a, 0)(^ Q x'>, 0) - 8 ua x'^a, 0)) J da^(y)F(y(a', r') - x(a, 0)) 

+ da^(y)5y x (a',r') £ — da^(z)F(z - y), (3.6) 

where P(X, fi, v) indicates a cyclic permutation of the indices A, /i, v inside the sum. The 
functional derivative of A with respect to x(er, 0) is then: 

J) A r 

fa „ (g 0) = (S m x"-(a,0)-d m x'"(a,0))Jda'"-(y)F(y(a',T')- x(a,0)) 
The quantity 5y x (a', r')/5x a (a, 0) depends on the type of surface that is considered. 
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Applying on both sides of Eq. ()3.7j) the operators d/dx'^(a,0) and d/dx 1 (a 1 0) in 
successive order, one obtains: 



d d 8A „ d 



da^(y)F(y(a',r')-x(a,0)) 



dx~f(a, 0) dx'P(a, 0) 5x a ((x, 0) dx^(a, 0) 

+/<^> (raw™?) S i/^'-'i' 



P(X,fi,v) 



(3. 



In order to satisfy the Bianchi identity (|2.11j) with representation f)2.15j) . the sum of the 
cyclic permutation of the left-hand side of Eq. ()3.8j) should vanish: 

V g 9 5A -0 (39) 

pfcp,*) 9 ^ 0)d*t>{a, 0)6x°(tr,0) ■ 1 • ' 

This involves through the right-hand side of Eq. (|3.8|) the sum Y^p(y,p, a ) af? aP^Jfs-- ^ or 
arbitrary contours, that quantity does not identically vanish. Therefore, the only solution 
to Eq. ()3.9j) is the constraint 

J2 T^x / da» u (y)F(y - x) = 0, x inside 5 or on C, (3.10) 

P(\,fj,,u) 

which is the defining equation of minimal surfaces. 

In order to check Eq. (|2.14jl . we limit ourselves to the case of minimal surfaces. In 
that case, the second term of the right-hand side of Eq. ()3.8)1 is null and the first term 
gives after contraction of the indices 7 and (3: 

_d d 5 A _ 1 x' a 

dxP dx'P 5x a ~ a (x' 2 ) 1 / 2 +(J ( a} - V- > 

On the other hand, regularizing the delta function of Eq. (|2.14j) according to the pre- 
scription 

one obtains, with the proper-time parametrization (x' 2 =constant): 

/l / 3 x la x llla x ,a x" 2 \ 

dy«5\y -x) = — [~ 2WW + ^ + ^ J + O(a). (3.13) 

Use of representation (J2.15j) and comparison of both sides of Eq. (J2.14j) (in its euclidean 
version) gives for finite a the identification (J2.16|) with C = 3/(4ir), which necessitates a 
quadratically vanishing g 2 N c with a when short-distance interactions are ignored. 
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Coming back to the first-order variation of the area [Eq. (|3.7|) ] it becomes in the case 
of a minimal surface: 



SA 

= (5, a x'»(a,0) - 8 va x'»(a,0)) I da^(y)F(y(a' ,t') - x(a,0)), 



5x a (a 1 0) 

(x ;2 x a — x'.xx >a ) 



x' 2 x 2 — (x'.x) 2 



(3.14) 



the second equation resulting in the limit a = 0; it manifestly satisfies the orthogonality 
condition 

X A 

in agreement with Eq. (|2.3jl . 

Some properties of minimal surfaces are presented in appendix [X] One main prop- 
erty that we shall use in the present paper concerns the commutativity property of two 
successive functional derivatives of the minimal area on its contour C, 

1 5 6 6 5 »)A(£7) = 0, (3.16) 



\5xP(a') 5x a (cr) 5x a (cr) 5xP(cr 
a feature that is reminiscent of a similar property of the Wilson loop average: 

4 The quark propagator in the external gluon field 

In dealing with the quarkonium bound state problem within the Wilson loop approach, 
one needs, in the path integral formalism, a representation of the quark propagator in the 
presence of an arbitrary external gluon field, satisfying the equation 

{i^f.d^x) ~ m ~ 91-A(x)^jS(x, x') = i5 4 (x — x). (4.1) 

The usual perturbative representation S(x, x') = i(i / -f.d — m — g / ~f.A)~ 1 5 4: (x — x') is not very 
convenient here since at each order of the perturbative expansion in the coupling constant 
the gauge covariance of the propagator is lost and the construction of the Wilson loop 
including contributions of fermion propagators becomes tricky. A representation that 
is well suited to exhibit the Wilson loop structure of the gauge invariant two-particle 
Green functions is the Feynman-Schwinger representation [22], which represents the quark 
propagator as a quantum mechanical path integral. This representation has however the 
main drawback that it dissolves the Dirac operator into the path integral and makes it 
difficult to obtain from it an equation that easily displays the properties of fermions. In 
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the present work we shall consider a representation that is based on an explicit use of 
the free Dirac propagator accompanied by a path-ordered phase factor. It will have the 
advantage of manifestly preserving the main properties of fermions. 

The building block of the representation is the gauge covariant composite object, 
denoted S(x,x'), made of a free fermion propagator Sq(x — x') (without color group 
content) multiplied by the path-ordered phase factor U(x, x') [Eq. ([2.1)1 ] taken along the 
straight segment xx'\ 

S(x,x') = So(x — x') U(x,x') . (4.2) 

The advantage of the straight segment over other types of line is that in the limit x' — > x 
U tends to unity in an anambiguous way. S satisfies the following equation with respect 
to x: 



(i T d {x) -m- gj.A(x)) S(x, x') = i5\x -x')+ z 7 a £ dX X {j^yj S (x - x'), (4.3) 

where the segment xx' has been parametrized with the parameter A as x(A) = Ax + (1 — 
X)x'. In the above equation, the point x' is held fixed; furthermore, the operator 5/5x a (X) 
does not act on the explicit boundary point x of the segment (corresponding to A = 1, cf. 
Eqs. (I2.2j) - (|2.3{l ). this contribution having already been cancelled by the gluon field term 
A. A similar equation also holds with respect to x', with x held fixed, with the Dirac and 
color group matrices acting from the right. 

The quantity —i(i~f.dr x \ — m — g^.A{x))8 i {x — x') is the inverse of the quark propagator 
S in the presence of the external gluon field A. Reversing Eq. 1)4.3)1 with respect to S^ 1 , 
one obtains an equation for S in terms of S: 

S(x, x') = S(x, x) - [ d 4 x"S(x, x")-i a f 1 dX X x 5 -ttt S(x", x'), (4.4) 



5x a (X) 

where the operator S/5x a (X) acts on the factor U of S, along the internal part of the 
segment x"x', with x' held fixed. Using the equation with x', or making in Eq. ()4.4)1 an 
integration by parts, one obtains another equivalent equation: 

S(x,x') = S(x,x')+ [ d 4 x" f 1 dX(l - X)S(x, x") - -. AN j a S(x", x'). (4.5) 
J Jo ox a (X) 

Equations f)4.4j) or (J4.5)) allow us to obtain the propagator S as an iteration series with 

respect to S, which contains the free fermion propagator, by maintaining at each order 

of the iteration its gauge covariance property. For instance, the expansion of Eq. ()4.4j) 

takes the form: 

S(x,x') = S(x,x') - ( d 4 y 1 S(x,y 1 ) 1 ai [' dX 1 X 1 —^ 7 —S{y 1 ,x') 



j jo 5x Ql (Ai; 

1 A A 

+ / 1 d\i x,dx 2 x 2 s{x, yih ai ( Sxai (Ai) s( yi , V2))i a2 ( 6xa2 (Aa) S(y 2 , x' 

+ ■••, (4.6) 
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the operator S/Sx(Xi) acting on the inner part of the segment yiUi+i of the phase factor 
U(yi, Vi+i)- A verification of Eq. ()4.1j) can be done through the above iteration series. The 
operator d/dx acts on three terms. First, acting on the free fermion propagator contained 
in S(x,x"), it gives, with the Dirac operator, a 5 4 (x — x") type term which then removes 
the integration over x" . Second, it acts on the boundary point x of U(x,x") and cancels 
the gluon field term A. Third, it acts on the inner part of the segment xx" of U. This 
term is then cancelled by the <5 4 -term coming from the next-order term of the iteration 
series under the action of the Dirac operator on the corresponding free propagator, and so 
forth. Equations (|4.4jl and (|4.5jl are generalizations of the representation used for heavy 
quarks starting from the static case [To] . 

The action of the operator 5/5x a (X) on U can be expressed in terms of an insertion 
of the field strength F [Eq. (|2.3J) ]. One can check with the first few terms of the series, 
using integrations by parts, that one can recover, in perturbation theory with respect to 
the coupling constant g, the conventional perturbative expansion of S in terms of g and 
A. 

5 The two-particle gauge invariant Green function 

The next step is to consider the quark- ant iquark gauge invariant Green function, for 
quarks q± and q 2 with different flavors and with masses m\ and m 2 : 

G(x 1 ,x 2 ;x' 1 ,x' 2 ) = (ip 2 (x 2 )U(x 2 , x 1 )'ip 1 (x 1 )tp l (x' 1 )U(x' 1 , x' 2 )tp 2 (x' 2 )) A:qim , (5.1) 

the averaging being defined in the path-integral formalism. Here, U(x 2 ,Xi) is the phase 
factor (j2.1j) along the straight segment x 2 x\ (and similarly for U(x[x 2 )). According to the 
conclusion reached in the final part of appendix El the dynamics of the system concerning 
its energy spectrum not containing string-like excitations can be probed by considering 
between the quark and the antiquark equal-time straight segments (in a given reference 
frame, e.g., the rest frame of the bound state); deviations of lines C X2X1 from the equal- 
time straight segment contribute only to the wave functional of the bound state and not 
to its energy. However, for covariance reasons, we first consider general straight segments 
x 2 X\ (not necessarily equal-time) and it is at a later stage that the equal-time limit will 
be taken. 

Integrating in the large N c limit with respect to the quark fields, one obtains: 

G(x 1 ,x 2 ;x' 1 ,x^) = -(ti c U(x 2 ,x 1 )S 1 (x 1 ,x' 1 )U(x' 1 ,x' 2 )S 2 (x' 2 ,x 2 )) A , (5.2) 

where Si and S 2 are the two quark propagators in the presence of the external gluon field 
and tr c designates the trace with respect to the color group. 
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The Green function G satisfies the following equation with respect to the Dirac oper- 
ator of particle 1 acting on x\. 



(i7-9(xi) -m 1 )G(x 1 ,x 2 ;x' 1 ,x' 2 ) = -i(tr c U(x 2 ,x 1 )5 A (x 1 - a/ 1 )Z7(a/ 1 , x' 2 )S 2 (x' 2 , x 2 )) A 

-27 Q (tr c f da(l-a) ^j^y \- Si(xi,x l )U{x' 1 ,x l 2 )S 2 (x 2 , x 2 )) A , (5.3) 
Jo ox a {a) 

where the segment x 2 Xi has been parametrized with the parameter a as x(a) = (1 — 
a)x\ + ax 2 ; furthermore, the operator 5/5x a does not act on the explicit boundary point 
Xi of the segment (cf. Eqs. ()2.2JI - (|2.3|) ). this contribution having already been cancelled 
by the contribution of the gluon field A coming from the quark propagator Si. A similar 
equation also holds with the Dirac operator of particle 2: 

G(xi,x 2 ;x'i,x 2 )(-ij.d( X2 ) ~ m 2 ) = -i(tr c U(x 2 , Xi)Si(xi, x[)U(x'i, x 2 )6*(x 2 -x 2 )) A 

+*(tr c f 1 daa S -^^-S 1 (xi,x' 1 )U(x' 1 ,x' 2 )S 2 (x' 2 ,x 2 )) A ^. (5.4) 
Jo ox p (a) 

Using for Si and S 2 representations (|4.4J) and (|4.5|1 . respectively, one obtains for G 
an expansion in a series of terms involving an increasing number of straight segments 
between xi and x[ on the one hand and between x 2 and x 2 on the other. With each 
segment is associated a path-ordered phase factor U ; the union of all such factors, together 
with U(x 2 ,Xi) and U(x' 1 ,x / 2 ), forms a Wilson loop along a skew polygon. We can then 
represent G in the following form: 

00 

G = ^ Gjj, (5-5) 

where Gij represents the contribution of the term of the series having (i — 1) points of 
integration between Xi and x[ (i segments) and (j — 1) points of integration between x 2 
and x 2 (j segments). We designate by Cij the contour associated with the term Gij. A 
typical configuration for the contour of G 4i3 is represented in Fig. |21 

Using for the averages of the Wilson loops appearing in the above series the represen- 
tation with minimal surfaces, and designating by A it j the minimal area associated with 
the contour Cjj, one obtains for the latter the following representation: 



G itj = (-l) l N c J d yi ■ ■ -d y^d z 1 ---d z i _ 1 S , 10 (a: 1 - Vih ai S 10 (y 1 - y 2 ) ■ ■ ■ 

X'f i - 1 Sio{y i -i - x'i)S 20 {x 2 - ^_i)7 ft - 1 5 20 (z i _i - z,-_ 2 ) • ■ ■ l Pl S 20 {z 1 - x 2 ) 
x£dn--- dn-i jT 1 dr[--- dr^_i(l - n) • • • (1 - n_x)(l - r[) ■ ■ ■ (1 - 7^) 

<0V1 m. (5.6) 



5x ai (ri) fe^-^Ti-i) fe^(r() Sxh-^r^i) 
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Here, the operators 5/5x(Tk) and 5/Sx(t^) act on the surface Aij through their action 
on the segments UkUk+i and z&i+i of the contour, respectively. (The parametrization of 
the segments of the line X\x' x is the opposite of that of Sec. EJ r now increases along 
the direction x\x' x to be in accordance with that adopted for the surfaces; cf. comment 
after Eq. (|3.1JI .) S w and 5*20 are the free quark propagators with masses mi and m 2 , 
respectively. 

Using in Eqs. ()5.3j) - (j5.4|) representations (|4.4jl - (|4.5jl for the quark propagators, one 
obtains: 



(i7A*i) - •m 1 )G(x 1 ,x 2 ; x[, x' 2 ) = -iS 4 (xi - x[) ^ G , 



.1 



00 r 1 5 



3=1 



00 

G(x 1 ,x 2 ;x' 1 ,x' 2 )(-i-f.d(x 2 ) - m 2 ) = +i5 4 (x 2 - x' 2 )^G ifi 

i=l 

"i . 5 



(5.7) 



r 1 



x(a)£xiX2 



(5. 



where G ,j and G^o are particular cases of representation (J5.6|) in which the particle 1 
or 2 propagators have been shrunk to a delta-function. These equations can also be 
obtained by making the Dirac operators act on Eq. ()5.5j) and using in the right-hand side 
representation (|5.bj) . 

The two equations (|5.3|) and (|5.4|) are independent, since the two Dirac operators 
concern independent particle variables. They are also compatible; this is evident from 
the very fact that G exists and is given by formula (|5.2|) . However, an independent check 
can also be done by making the Dirac operator of particle 2 act on Eq. ()5.3jl and the 
Dirac operator of particle 1 act on Eq. (|5.4j) . using in the right-hand sides the properties 
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of the propagators and of the phase factors U and subtracting from each other the two 
resulting equations. The result is zero, due mainly to the facts that the final expressions 
involve functional derivatives 5/Sx(a) and 5/5x(a') of the phase factor U(x2,Xi) and 
these commute. A similar verification can also be done with Eqs. (I5.7jl - (|5.8j) . where 
now the representation of the Wilson loop averages by minimal surfaces has been used. 
One has to use in the right-hand sides properties of the terms G^ and cancellations 
between contributions of successive GijS along the quark and antiquark lines. One ends 
up with expressions involving functional derivatives 5/5x(a) and 5/5x(a') of minimal 
surfaces along the segment x±X2- The compatibility of the two equations is then due 
to the fact that two successive functional derivatives of a minimal area localized on its 
contour commute [Eq. ()3.16|) and appendix IX] . 

Eqs. (I5.7|) - (j5.8|) are not closed integro-differential equations for G, for once the action 
of the functional derivatives 5/5x on the various minimal surfaces has been evaluated one 
does not obtain back G on the right-hand sides. This feature makes difficult the search 
for a bound state equation in compact form. In this respect, if G has a bound state pole 
in momentum space, then the right-hand sides of Eqs. ()5.7|) - (j5.8|) should also have the 
same pole; this is possible only if the actions of the functional derivatives 6/5x on the 
partial ingredients Gij of G yield among other terms common factors that allow coherent 
summations of the Gijs to produce again a pole term; otherwise, each Gij, containing a 
finite number of free quark propagators, cannot produce alone such a pole. In x-space, 
the selection of a bound state is made by taking a large time separation between the pairs 
of points (x\, £2) on the one hand and (x[, x' 2 ) on the other |28j . 

The independence of the two equations ()5.7p - (j5.8j) means also that they might allow 
the elimination of the relative time variable of the two particles prior to any resolution of 
an eigenvalue equation and the reduction of the internal dynamics to a three-dimensional 
space, a feature which was outlined at the beginning of this section according to the re- 
sults obtained at the end of appendix |XJ However, such a reduction does not seem easily 
manageable on the general forms of Eqs. (|5.7|) - (J5.8J) . Furthermore, arbitrary approxi- 
mations made in the right-hand sides of those equations may destroy their compatibility 
property. This is why in the following we shall directly study the equal-time limit of the 
system by considering the "sum" of the two equations and then, at a later stage, shall 
indicate how to determine its relative time evolution law by considering the "difference" 
of the two equations. 
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6 Bound state equation 



In the large N c limit, the mesonic sector of QCD is composed of one-particle states, which 
are bound states of a quark- ant iquark pair and of gluons fHl E3J Taking in the Green 
function ()5.1|) a large time separation between the pairs of points (x%, x 2 ) and (x\, x 2 ) and 
using the completeness relation one finds: 

G(x 1 ,x 2 ;x' 1 ,x' 2 ) =Y J ^n(x 1 ,x 2 )^n(x' 2 ,x' 1 ), (6.1) 

n 

where $ n is the wave functional of the bound state labelled with the collective quantum 
numbers n: 

i^n, ai ,a 2 (x u x 2 ) = < 0\tp 2a2 (x 2 )U(x 2 ,x 1 )ip hai (x 1 )\n >, 
-^,^04, xi) = < n\^ la , i (x' 1 )U(x' 1 ,x 2 )ij 2 ^ 2 (x' 2 )\0 > . (6.2) 

Since the lines C X2Xl and C x ^ x i are rigid straight segments completely determined by 
their end points, one can consider the wave functionals as functions of the end point 
coordinates. Introducing total and relative coordinates and momenta, 

X = ^(x 1 +x 2 ), X = X 2 -Xi, P = Pl +p 2 , p=-(p 2 -p 1 ) ) 
d_ 

OXa 

and considering a bound state with total momentum P, one has: 

$(xi,x 2 ) = e~ iRX ^P,x). (6.4) 

In the large separation time limit, the right-hand side of Eq. (J6.1j) displays a series 
of oscillating functions in the separation time variable. By appropriate projections and 
integrations one can select in this series the bound state that will survive in the large 
time limit [2E] • It should be emphasized that in the above limit only terms that factorize 
in G into expressions depending on the line 2:1X2 and expressions depending on the line 
x[x 2 could survive to the selection operation of the bound state. Terms that still contain 
expressions joining line x[x 2 to line X\X 2 would not contribute to the previous operation 
and could be discarded. 

It is convenient to consider the total momentum P of the selected bound state as 
a reference timelike vector and define transverse and longitudinal parts of vectors with 
respect to it: 

€ = %- q ~§iP^ qi = (q.P)P„ h = ^=, qL = q-P, 



Pa, l M = i-^i, a =1,2, (6.3) 



Pl = VP 2 , q L =q , q T2 = -q 2 , V ' -x T2 = r. (6.5) 



2 ,/_^T2 
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These decompositions are manifestly covariant. To further simplify the notation we adopt 
for the Dirac matrices the following convention: they will be written on the left of the 
spinor functions with labels 1 or 2 indicating on which particle indices they act, the 
particle 2 matrices (the antiquark at x 2 ) acting actually from the right; thus: 

l2^il2uG = G ai p 2!(x 1 

Similar definitions also hold when G is replaced by $. Notice that products of 72 matrices 
act from the right in the reverse order. 

We next introduce the free Dirac hamiltonians of the two particles: 

ho = mi7ii - JiliT-pI, ho = -m 2 72L - l2Lll-pl- (6.7) 

Going back to Eqs. ()5.7j) - (j5.8j) . multiplying the first by 71L, the second by — 72L, adding 
and subtracting the two equations one finds: 

00 00 
(Pil + P2l) ~ (ho + h 20 ) G = -iS A (x 1 - x' 1 )7i L G o,j ~ ^ A (x 2 - x' 2 )j 2L G i,o 

j=l i=X 



+i 



7lL7l ° So 1 Ml - a) 6^) + l2Ll " Jo daa jJ^ I £ G 



1,3 



, (6.* 



(Pil ~ P2l) - (ho - ho) G = -i5 4 (x 1 - x[)-f 1L ^ G ,j + iS A (x 2 - x' 2 )^ 2L G ifi 

j=x i=i 



+i 



. (6.9) 



Equation ()6.9|1 mainly determines the relative time evolution of the Green function, while 
Eq. (|6.8|) mainly determines the dynamical properties of the two-particle system. Since 
that equation does not involve the relative energy operator (p 2L — p iL ) in its left-hand 
side, one is entitled to take in it the equal-time limit xl = 0; the equation becomes in 
that case three-dimensional with respect to the transverse relative coordinates x T . After 
determining, within a given approximation, the dynamical properties of the system from 
Eq. ()6.8|) in the limit xl = 0, one can go back to Eq. ()6.9|) and determine, within the same 
approximation, its relative time evolution law. We shall henceforth follow this method of 
approach. 

In the large separation time limit between the pairs of points (x\,x 2 ) and (x[,x' 2 ) 
the delta- functions that are present in the right hand-sides of Eqs. ()6.8|) - (|6.9J1 do not 
contribute and can be ignored. In order to extract from the right-hand side of Eq. (|6.8|) a 
bound state wave function $, it is necessary that the actions of the functional derivatives 
S/Sx(a) on the various parts Gij yield among other terms common factors that factorize 
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the Green function G again. The functional derivative S/5x(a) acts on a given Gij [Eq. 
(|5.6|) ] through the exponential factor containing the minimal area term and yields a term 
proportional to 5Aij/5x(a). This term may itself be acted on by the other functional 
derivatives existing in the definition of Gij. One thus ends up with two types of term. The 
first contains the first-order derivative 5Aij/5x(a), which factorizes the other derivatives 
on the left, and the second contains terms involving higher-order derivatives of A^j, one 
of the derivatives acting along the straight segment X\Xi. The dominant part in the large- 
distance limit comes from the first type of term; higher-order derivatives of the minimal 
area tend to weaken the large-distance behavior (cf. appendix [X] for the second-order 
derivative). In the following we shall mainly concentrate on the contribution of the first 
type of term; terms containing second-order derivatives will be considered in Sec. [71 in 
connection with the self-energy parts of the quark propagators. 

The term that is retained with Gij is an integral over a of a function proportional 
to 5Ai > j/Sx(a). According to Eq. (j3.14|) . the latter depends only on the local properties 
of the surface at the point x(a, 0) lying on the straight segment x±X2, namely upon the 
derivatives x'(a,0) and x(a, 0). Because the line C X2Xl is a straight segment, one has 
x'(a, 0) = x, independent of the form of the surface. In the equal-time limit (xl = 0) 
the previous relation becomes x'(a, 0) = x T . The derivative x(a, 0) depends, however, on 
the form of the minimal surface in the vicinity of the straight segment. Using methods 
of analysis similar to those used at the end of appendix El one can show that in the large 
separation time limit x(a, 0) tends to a linear function of sigma. In that case, one can 
parametrize it as x(a, 0) = (1 — cr)xi + ax2, where x\ and £2 are the slopes at the points 
X\ and x 2 , respectively, ii = ±(0,0), ± 2 = ±(1,0). The latter make still reference to 
the other end points of the corresponding segments; for the case of the simplest contour 
Ci 5 i, these are x[ and x' 2 . In order to remove any explicit reference to the points of the 
remote past an operator representation of the slopes, depending only on points x\ and X2 
becomes necessary. 

To find such a representation, we consider the simplest contour C^i corresponding to 
Gi t \. Here, one has X\ = (x[ —Xi), ±2 = (x' 2 —x 2 ), x a i < with our parametrization and 
because of the facts that x a i — > +00 and x' aL — > —00 (a = 1,2). Considering in general 
the cases xl and x' L finite in the above limits, one can set, modulo negligible terms, 
Xil — X 2L- (This is exact in the equal-time cases xl = and x' L = 0.) Furthermore, 
a close examination of the integrals of the term 6Ai s i/6x(a) shows that the x a terms 
(a = 1,2) appear after integration in the forms x"/\x a L\ and through their orthogonal 
components to x±X2- An operator representation of the latter can be found by making 
the quark momentum operators p afl act on the term G\^\. They generally yield three 
different contributions. The first comes from their action on the corresponding free quark 
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propagator, the second from the segment X1X2 and the third from the segment (x' a — 
x a ). Since the terms x"/\x a L\ already appear in expressions that are proportional to the 
string tension a, one can use for the latter terms, as a first approximation, free theory 
expressions; in this case, it is sufficient to retain the contributions coming from the free 
quark propagators. The latter, in the large time limit yield with massive quarks the 
following dominant behavior: 

p af ,S a o{x a - x' a ) ~ m a f tt X *=£ s S aa (x a - x' a ), (x a - x' a ) L -> +00, a = 1, 2, (6.10) 

(x a - x' a ) 2 



from which one deduces: 



1,2, (6.11) 



\&ol\ \PaL\ 

where the operators p a L are the free theory expressions of the individual energies: 



PaL = h a0 , \p aL \ = yjml -pi 2 = E a (pT), a = 1,2, (6.12) 

h a0 being defined in Eqs. (16. 7}) . 

This approximation will retain all terms of order up to 0(l/c 2 ) in the nonrelativistic 
limit; terms that have been neglected, if they are nonzero, have contributions in the 
nonrelativistic limit starting at order 0(l/c 4 ). 

We next generalize the above evaluation to the higher-order contours Cij appearing 
in Gij. Here, we make the assumption that in the large separation time limit between the 
pairs of points (x\, X2) and (x[, x' 2 ) the derivatives SAij/ 5x(a) can be expanded around the 
driving term 8A\ : i/8x(a) of the lowest order surface. Neglecting the higher-order terms 
of these expansions, one ends up with the common operators 8Ai^/5x(a) to all factors 
Gij, involving the segment X\X 2 and the momentum operators p afl , a = 1,2. Those terms 
can then be factorized in front of G and interpreted as potentials. 

The bound state equation obtained from Eq. ()6.8|) in the equal-time limit xl = is 
then: 

Pl - {h w + h 20 ) - IilI^A^ - ^l^A 2 ^(P l , x T ) = 0, (6.13) 
where ip is the wave function in the equal-time limit, 

ij(P L , x T ) = <j>(P L , x L = 0, x T ), (6.14) 

and the potentials are defined through the equations (in Minkowski space) 

^/Vd-O^j, A^oJ^a'o'^. (6.15) 

They can be calculated either by using the minkowskian version of Eq. (jSH)) and then 
conditions ()6.11|) - (j6.12|) . or by first writing Eqs. ()6.8j) and ()6.13|) in euclidean space, using 



21 



Eq. (|3.14|) and then passing to Minkowski space. (In euclidean space the right-hand sides 
of Eqs. f)6.15j) contain an additional (— i) factor.) Since 5Ai t i/5x(o J ) is orthogonal to x 
[Eq. ()3.15j) ]. the resulting vectors will satisfy this property. We define transverse vectors 
with respect to x with a superscript "i" : 



Qu r,X.Q. 

X 1 



(6.16) 



However, x itself is orthogonal to the total momentum P in the equal-time limit (xl = 0) 
and reduces to x T . The part of the three-dimensional relative momentum p T that is also 
orthogonal to x T will be denoted p Tt : 



Tt 



T 



T T T 



(6.17) 



This vector enters in the definition of the relative orbital angular momentum. Defining 
the corresponding Pauli-Lubanski vector Wl (L referring here to the orbital angular 
momentum) as 



one has the relations 



wl 



^0123 — +1; 



-P 2 [x T2 p T2 ~ x Ta x T Pp T a p T p - 1ix T .p T ) = -P 2 x T2 p Tt ' 
The expression of A 1 ^ l is: 



W 2 



A 



-a\ -x 



„T2 



E\E 2 



E 1 + E 2 



9^(p 2L ) - e(pii)) 
x T2 P 2 E\E 



E 1 E 2 



p2 



,x T2 P 2 



(6.18) 



L 2 . (6.19) 



(E 1 



E 2 f y 2W 2 



■e x e 2 - r 



2 Tt 
P 




arcsm 



(E 2 + E\ 
-x T2 P 2 



\{,e{p 2L ) - e(p 1L )) 
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E 



-w 2 L 



-x 



T2 



p2 



Tt 

pI 



- Ei + E 2 - 



-Wl 
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-w 2 



-x T2 P 2 E 2 



-Wl 



\ -x T2 P 2 E 2 




i \\ EiE 2 Tt 
e{pi L )) - ■ - +p L 
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Ei 



-x 



T2 



-wl 



-x T2 P 2 E 2 



(6.20) 



*P*E%^ E 1 + E 2 \ 

Here, e(pn) and e{p 2 L) are the energy sign operators of the free quark and the antiquark, 
respectively: 

e( PaL ) = ^, a = 1,2, (6.21) 



E n 
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h a o and E a being denned in Eqs. (J6.7|) and (|6.12|) . The expression of A 2fI is obtained from 
that of Ai^ by an interchange in the latter of the indices 1 and 2 and a change of sign of 
p Tt . In particular, the longitunal parts of the potentials A\ and A 2 add up in Eq. ()6.13j) . 
One has for their sum the expression: 

E\E 2 [(Ei E 2 , 




+ arcsm 




-wi 



-x 



T2 



p2, 



+ (e(pii) - c(P2l)) 



-Wl 



-wi 



-x T2 P 2 E 2 \ -x T2 P 2 E\ 

(6.22) 



For sectors of quantum numbers where W 2 = 0, the expressions of the potentials 
become: 



A 1L + A 2L = l -{e{ PlL ) + e(p 2L ))(nf^2~, (6.23) 



< = -Ek^ {Ei+E2) 'l Ei >>^ 

K = + ^e- 2 (1 {Ei + E2) " l E2 ) p >^-^ 2 - (6 - 24) 

All expressions of the potentials have been written as classical functions of their argu- 
ments, without taking into account ordering problems. These necessitate a detailed study 
which will not be done here. We simply outline some general features that may be useful 
for the resolution of the wave equation, (i) Many ordering problems that concern linear 
momentum operators do not affect the energy eigenvalues and rather concern the defi- 
nition of the kernel of the scalar product of the wave functions; one can pass from one 
definition to the other by appropriate changes of function, (ii) The square-root and arcsin 
functions which involve the variables x T2 and 1/E 2 (a = 1,2) could be treated in first ap- 
proximation by replacing in E a the radial momentum operator squared by its mean value 
in the bound state, or, if the resolution is done in momentum space by replacing x T2 by 
its mean value, (iii) A close study of the chiral properties of the wave equation suggests 
us to further adopt the following rules: the doubly transverse momentum operator p^, 
maintaining its definition of Eq. (|6.17|) . and the energy sign operators e(p a L) ( a = 1,2) 
should be placed on the utmost left. 

Equation (|6.13j) . together with expressions (J6.20J) . is very similar to an equation pro- 
posed by Olsson et al. on the basis of a model where quarks are attached at the ends of 
a straight string or a color flux tube [33J EH] ; the difference mainly concerns the energy 
sign operators; in Ref. the equation which was proposed is the Salpeter equation 



23 



|3Uj . in which generally the potential is proportional to a global energy projecter; here, 
the energy sign operators, though they could be expressed through energy projectors, do 
not match exactly the projector of the Salpeter equation; the doubly transverse parts 
of the vector potentials do even not have energy sign operators. Apart from this slight 
difference, however, the physical significance of Eq. (jOHj) is the same as that of Ref. [31] . 
The vector potentials A afl (a = 1,2) can be interpreted as representing contributions of 
the energy- momentum vector of the straight color flux tube with variable length r; its 
energy is represented by the sum A 1L + A 2L [Eq- (|6.22j) ]. while its angular momentum 
contributes through the doubly transverse components A™. 

The norm of the wave function ip can be obtained (after a few approximations) from 
Eq. ()6.8|) . The details are presented in appendix [Bj The result is: 



where h a and E aj a = 1,2, are the free Dirac hamiltonian and energy of each particle, 
including now the self-energy contributions [cf. Sec. [7] and Eqs. (I7.9j) - (j7.11|) ]. 

Finally, let us mention that Eq. (J6.9j) can be used to determine the relative time evolu- 
tion of the wave function $ [Eq. (J6.2j) ]. Using in the right-hand side of Eq. (J6.9j) the same 
types of approximation as in Eq. ()6.8|) . taking the large separation time limit between the 
pairs of points (x\, x 2 ) and (x[, x' 2 ) and passing to the bound state wave function, one can 
integrate Eq. (|6.9|) for the latter obtaining the relative time dependence in the form of 
an ordered exponential function involving the various operators and potentials appearing 
in the equation, the initial condition being given by the function iP(Pl,x t ) [Eq. (j6.14j) . 
We shall not need, however, that expression in the present work. 

7 Quark self-energy 

In general, there are functional relations between kernels of two-particle Green func- 
tions and self-energies of the constituent particles j3JJ|. In conventional Green function 
equations (defined without phase factors), one easily factorizes the self-energy contribu- 
tions and incorporates them into the external propagator contributions. Furthermore, 
self-energy parts can themselves be evaluated by means of Schwinger-Dyson equations. 
In the case of gauge invariant Green functions, however, the factorization of self-energy 
contributions becomes a hard task, since the path-ordered phase factors and the result- 
ing Wilson loops maintain all interacting pieces into global entities. The presence of 
self-energies is, however, necessary for a consistent study of the properties of the sys- 
tem under given symmetries, transcribed usually into Ward-Takahashi identities. In the 
present case, the symmetry that is of interest is chiral symmetry. Also, when the dynam- 




(6.25) 
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ical equations of confining interactions are expressed in momentum space, the presence 
of self-energy contributions becomes necessary to remove the infrared divergences from 
observable quantities [37] . 

When the interaction kernel of the bound state equation is represented by the me- 
diation of an effective propagator, one expects to obtain the self-energy contribution by 
contracting the lines of the outgoing particles (assumed to be of the same type) through 
a single particle propagator and forming a loop with the kernel. In the present case, the 
interaction part of the bound state equation (J6.13|) can be visualized by multiplying it 
back by the factor 7il72L- It has three different tensor parts: the first corresponds to ver- 
tices with the matrices 7il72l! the second and third to vertices with matrices 7p.pf*72L 
and 7ii7| , *-P2 > respectively. When these terms are incorporated in a two-point loop and 
integrated in momentum space, the non-invariant pieces under spatial rotations disappear 
and what remains is simply the part of the interaction with the 71L72L matrices, inside 
which also the angular momentum operator has disappeared (cf. Eq. (I6.23j) ). This cor- 
responds to the situation where the interaction kernel is generated by the mediation of a 
Coulomb-gluon propagator, proportional in x-space to 5(xl)V— x T2 . 

The above result also coincides in form with that obtained in two-dimensional QCD. 
Here, in the large N c limit, one has two different but equivalent ways of obtaining the 
bound state equation: either by working with the Schwinger-Dyson approach in the 
axial gauge, or by working with the Wilson loops with minimal surfaces jT^J. (It is 
even easier to work in the light cone gauge |3Z|, but the latter is less useful for our 
purpose here.) The Bethe-Salpeter equation with the one-gluon exchange kernel and 
the corresponding self-energies becomes an exact bound state equation and because of 
the instantaneity of the propagator, it yields the Salpeter equation [20]. The procedure 
developed in the present paper is also applicable to two-dimensional QCD. It is sufficient 
to remove from the results that were obtained the rotational motion part. In this case, 
the interaction potential reduces to the expression of Eq. (jfi.23)) which is nothing but 
Salpeter's kernel. What misses however in the bound state equation is the corresponding 
self-energy contributions. We shall now show that Eq. (J6.8J) . which is at the origin of the 
bound state equation, also yields the quark self-energies. 

Since the formal expressions of the self-energies are the same in both two and four 
dimensions in x-space, we shall directly work in two dimensions, by freezing the transverse 
variables with respect to x and considering the simplified case of a plane. A self-energy 
contribution for particle 1 will be recognized as depending only on the variables of that 
particle, not making reference to the variables of particle 2, except for the directions of the 
total momentum P of the bound state, which is chosen as the reference timelike direction, 
and of the relative coordinate (x 2 — Xi). The simplest contribution to the self-energy of 
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particle 1 comes from the term G 2 ,\ in Eq. ()6.8|) . a typical contour of which is shown in 
Fig. Ela. 





(a) 



(b) 



Figure 3: (a) A typical contour 6*2,1 associated with the term G2,i of the two-particle 
propagator, (b) A configuration where the segments X1X2 and y~\x' x intersect. 

The term that is relevant here is that in which the operator S/Sx(a) (multiplying the 
71 matrices) acts on the term 5A2,i/5x(t) of G 2j i [Eq. ([5.6)1 ] . One thus obtains the second- 
order functional derivative 5 2 A 2) \j 1 5x^{a)5x a {r) , where r parametrizes the segment y\x\ 
and a the segment X\X 2 . The expression of the second-order functional derivative of 
the minimal area has been given in Eq. (jA.l|) . where now x[o\ 0) and x(a, 0) have to be 
replaced by x(a) and x(r), respectively. Since x(a) and x(r) belong to different segments, 
the terms proportional to the explicit delta-functions can be dropped; furthermore, since 
the orthogonal variables to the surface have been frozen, also the last term of that equation 
can be ignored; it is only the second term of the right-hand side of the equation that may 
contribute: 



where x (r) andx(a) are parametrized as x{r) = yi+r(x' 1 —yi) andx(cx) = xi + a(x 2 — xi). 
In the limit when the regulator a vanishes, the function F, Eq. ()3.3|) . tends to a two- 
dimensional delta- function [Eq. (|3.4|) ]: therefore, the above second-order derivative is 
nonvanishing only when the two segments are intersecting (Fig. [31b) ; we then have to 
integrate with respect to a and r with the weight factors (1 — cr) and (1— r) [Eqs. f|5.6|) and 
(|6.8|) ]. The calculation can be done by replacing F by a two-dimensional delta-function 
involving a and r and the arguments of which can be obtained by first writing the explicit 
expression of (x(cr) — x(t)) 2 with two-dimensional components x° and x 3 (in the bound 




x 



(5 Xa (x 2 - x x ) v - 5 ua (x 2 - a; 1 ) A )((J AAt (x , 1 - y x ) v - 5 u ^x\ - Vl ) x ), (7.1) 
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state rest frame) for instance: 

F(x(a)-x(T)) = 5(x°(a)-x°(T))5(x 3 (a)-x 3 (T)) 
= S(z° + ax° - ry°) 5(z 3 + ax 3 - ry 3 ) 



2x 3 y 3 



8 8 

crar - ry") — — \z 6 + crrr 
da or 



ry 



(7.2) 



where we have defined x 



x 2 - Xi, z 



xi - yi, y 



yi. Actually not the full 



expression of the above integral is needed, but rather that part which depends only on 
the points X\ and y\, which correspond to the integration endpoints a = and r = 0. 
Integrating that expression by parts and retaining only the latter endpoints, one obtains: 



£ dadr (1 — cr)(l — t)f{x{<j) - x(r 



1 



2x 3 y 



-5(z°)\z 3 \ + ... 



(7.3) 



where the dots stand for the remaining terms that do not contribute to the self-energy. 
The contribution of the integral of the second-order derivative of the area, Eq. (|7.1|) . is 
then: 



/ dadr (1 - a)(l - t) 
Jo 



8A 2 a 



5 x>* (a) Sx a (r) 



-<VxU5(z°)|* 3 



S^LS aL S((x 1L - yi L )y-(xi - yT) 2 



(7.4) 



where we have dropped additive terms not depending only on z and have done the tensor 
calculation in two dimensions and restored at the end the covariant expression. 

The presence of the delta-function along the temporal direction means that the seg- 
ments X\X 2 and y±Xi are parallel, or more generally lie in an orthogonal plane to the time 
direction (Fig. 




Figure 4: Configuration of the contour C24 corresponding to the self-energy contribution. 



In the limit (xil — x' 1L ) —> 00, the areas A 2 ,i (x' 2 X2Xiyix' 1 x' 2 ) and A 1 i(x 2 X2yix' 1 x 2 ) 
become almost equal and one may replace in Cr2,i the former by the latter. The expression 
(|7.4j) . combined with the multiplicative free quark propagator Sio(xi?/i), then factorizes 



27 



G*i i and plays the role of a self-energy correction. In order to complete the derivation, 
one must repeat the same calculations with the higher-order terms Gij (i > 3, j > 1), 
where in the second-order derivatives of the type of 1)7.1)1 . the functional derivative 5/Sx(t) 
corresponds to 5/5x(ti) and acts on the second segment 7/17/2 of the corresponding contour 
(Fig. EJ). One thus finds the first-order self-energy correction (in the string tension 
a) that factorizes the bound state wave function in Eq. ([6.13)1 : it corresponds to the 
effective exchange of a Coulomb-gluon, with propagator proportional in momentum space 
to 5, L 5 aL l/{{-q T ff. 

The higher-order self-energy corrections are obtained by extending the above procedure 
to terms that contain higher numbers of second-order derivatives of the areas Aij. The 
next correction is provided by the terms G{j with % > 4 and containing two second-order 
derivatives of the area. Let us consider for definiteness the term G4.1 corresponding to 
the contour C^i (Fig. 0a). 




Figure 5: (a) A configuration of the contour C^i corresponding to the second-order self- 
energy correction, (b) A configuration corresponding to overlapping integrals. 



The self-energy contribution comes from the term containing the product sx(a)5x(T 3 ) 
5x(t^)5x{t 2 ) i where S/5x(n) acts on the segment t/it/ 2 , 5/5x(t 2 ) on y 2 y 3 and 5/5x(t 3 ) on 
7/32/1. The resulting delta-functions imply intersection of the segments X\X 2 and 7/3X1 
and parallelism of the adjacent segments 7/17/2 and 7/22/3. The two delta-functions can 
be integrated as before in an independent way. The result is proportional to the prod- 
uct 5((x 1L - yzL)\J-{xT - t/J) 2 S((y 1L - y2L)\J-(yT - yi) 2 - Notice that in the case of 
overlapping integrals, corresponding to the situation where two derivatives act on the 
segments Xix 2 and 7/27/3 on the one hand and on the segments 7/17/2 and 7/3X1 on the other, 
the delta- functions would imply intersections of the segments of each of the above pairs 
(Fig. 0b); integrating the delta-function of the second pair as before yields the delta- 
function 5(t/2l ~ 2/3l); the implementation of which prevents the segments XiX 2 and 7/27/3 
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from intersecting and the corresponding integral vanishes. Also notice that in the two- 
dimensional limit that we are considering here (freezing of orthogonal deformations of the 
areas) higher-order derivatives of the areas vanish, for they would involve derivatives of 
the function F outside its support. 

The generalization to higher orders is now straightforward. Any high-order contribu- 
tion will contain in the term Gij a product of second-order derivatives of the area A it j, 
in which one of them is Sx ^ Sx ^ Tk ^ (k < i — 1), where 5/5x(rk) acts on the segment ykUk+i] 
the others are nonoverlapping (disjoint or nested) and act on segments lying between x\ 
and yk+i, being thus nested within xii/k+i. The remaining part of the calculation repeats 
the steps used for the first-order correction. 

One thus generates the whole series of self-energy corrections that can be summed 
into the Schwinger-Dyson integral equation, corresponding to a Coulomb instantaneous 



kernel -a5 lxL 5 1/L 5(x L - y L )\J-{x T - y T ) 2 . 

Defining the momentum space Fourier transform of the function \J — x T2 through an- 
alytic continuation of the power parameter , one finds: 



/rfW^ T v^=-7-^. (7.5) 

Denoting by S(pL,p T , m) the quark propagator in momentum space with free mass m tied 
to the bound state of momentum P and defining the self-energy contribution through the 
equation 

iS~ 1 (pl,P , rri) = j.p — m — S(p T , m), (7.6) 
one has the Schwinger-Dyson equation: 

f d 4 k 1 
E(/,m) = 8trx 7l J j^S(k L , k T , m) _ fcr)2)a 7L. (7-7) 

Because of the instantaneity of the kernel, the self-energy E actually depends only on 
the three-dimensional transverse momentum p T . The tensor decomposition of E can be 
written in the form 

£(p T , m) = 7 T .p T A(-p T2 , m) + B(-p T2 , m). (7.8) 

The self-energies of the quark and of the antiquark have now to be incorporated in the 
Dirac energy operators in Eq. (|6.13J) and in the definitions of the energy sign operators 
f[6~2T|). Defining 

hiip[,mx) = h 10 +j 1L E(pJ,m 1 ) = 7u(mi -n/T-Pi + S(pf,mi)) 

= 7l£ ((m 1 + SO- 7^(1-^)), 
h 2 (p2, m 2 ) = h 20 - 72lE(-^, m 2 ) = -7 2 L(m 2 + 7J .p 2 + , m 2 )) 
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-l2 L ((m 2 + B 2 ) + 7 2 T -^(1 - A 2 )) , (7.9) 



E 2 {-p T 2 \m 2 ) = y/h% = V(m 2 + J B 2 ) 2 -pP(l-A 2 ) 2 , (7.10) 

K 

Ea 



e( PaL ) = ^, a = 1,2, (7.11) 



equation ()6.13|) becomes 

Pl ~ {hi + h 2 ) - 7 il7i ^ - 72L72^]^(^L, x T ) = 0, (7.12) 

where the potentials A are given by expressions (J6.20J) in which E a and e a (a = 1, 2) are 
replaced by expressions (|7.10|) and (j7.11j) . respectively; furthermore, the self-energy func- 
tions should also be incorporated in the appearances of the doubly transverse momentum 
p Tt and the orbital angular momentum Wl. Since we have replaced p\ l and p\ l in favor 
oip Tt , the corresponding substitutions may not seem straightforward. However, the dif- 
ficulty is circumvented easily. Each energy factor E a (a = 1, 2) that appears explicitly in 
the potentials is reminiscent of a term of the type / E a ; p\ l undergoes the substitution 
pj* — > (1 — Ac^Pq 1 , where A a (a = 1,2) is the self-energy function, Eqs. (J7.9j) . It is then 
sufficient to replace each E a in the potentials by E a /(1 — A a ), without modifying the 
momentum and orbital angular momentum operators. 

The quark propagator is not a gauge invariant quantity and therefore could not lead 
to observable effects. This feature manifests itself through the self-energy equation (|7.7jl . 
which, due to the infrared singularity of the integrand, provides an infrared divergent self- 
energy. On the other hand, physical observables should be free of infrared singularities. 
The treatment of infrared singularities in momentum space depends on the method of 
regularization that is adopted. Dimensional regularization, which has the advantage of 
preserving symmetry properties, and which we adopt throughout this work, gives zero for 
the abovementioned singularities, but has the disadvantage of hiding the distinction be- 
tween unobservable and observable quantities. This is why it is necessary to check at every 
stage of calculations whether physical quantities are indeed free of infrared singularities 
by also using an explicit infrared cutoff method. 

In the present approach, we have to check that the wave equation (|7.12|) . which de- 
scribes the physical properties of the bound states, is free of infrared singularities. To 
do this, we have to isolate the singularities coming from the self-energies. The properties 
of the latter will be studied in Sec. |H1 here, we mention the most relevant ones for our 
purpose. Designating by V the linear confining potential, Eq. (|8.3|L it is found that the 
singularities of the self-energy parts are contained in the energy factors E a [Eqs. (J7.10|) ]. 
while the ratios (1 — A a )/(m a + B a ), a = 1,2, are free of singularities; from the integral 
equations satisfied by E a , Eq. it is seen that their singularity is represented by 
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the three-dimensional integral in momentum space of — V/2. [Section |H] deals with the 
equal- mass case and mainly with the massless limit; however, the infrared singularity 
properties are not affected by the values of the quark masses and could be abstracted 
from the equations of that section.] Therefore, the combinations (E a + V/2) are free of 
infrared singularities. Coming back to the potentials A a ^ [Eqs. (|fi.2()j) ] their large- distance 
behavior in rr-space can be studied by expanding the various functions contained in their 
expressions in terms of (1 — A a ) 2 ~L 2 / (E 2 x. 2 ) , a = 1,2 (in the cm. frame). Factorizing 
the function or, the expansions are of the type ar(l + 0((1 — A a ) 2 L 2 / (E 2 x. 2 ))); therefore 
the higher-order terms are nondominant at large distances and could not lead to infrared 
singularities. For the present study it is sufficient to keep the leading terms, which are 
represented by the expressions ()6.23j) - ()6.24|) . 

Considering first the contribution of (An + A21) and replacing the energy sign op- 
erators by their expressions (j7.11|) . we immediately find that the energy factors and the 
potential appear with the combinations (E a + V/2), which ensure the infrared finiteness 
of the result. Next, we consider the contributions of the spacelike potentials A^ [Eqs. 
(j6.24j) . where each E a should be replaced by E a /(1 — A a )]. The doubly transverse momen- 
tum operator p Tt is equal in the cm. frame to p* = — x L. The effective potential that 
matters is proportional to arp* = — x L. The orbital angular momentum operator L 
acts on spherical harmonics and does not modify the infrared properties of wave functions. 
The term x/r has as Fourier transform a function proportional to p/p 4 , which does not 
lead to infrared singularities, after the angular integrations in convolution integrals are 
done. Therefore, the wave equation ()7.12)) is free of infrared singularities. Explicit cases 
of the above cancellations can also be found in the equations presented in Sec. |H1 

8 Chiral symmetry breaking 

The presence of the self-energy contributions in the bound state equation ()7.12|) allows 
us to study the possibility of the spontaneous breakdown of chiral symmetry. This is 
intimately related to the existence of nonperturbative solutions to the Schwinger-Dyson 
equation in the chiral limit when the quark mass m tends to zero |3U| EE] • 

Using decomposition ()7.8|) and considering the bound state rest frame, the integral in 
Eq. ()7.7|) can first be integrated with respect to the energy variable ko giving rise to two 
coupled equations: 



B(p) 



(k 2 (l- A(k)) 2 + (m + B(k)) 2 ) 
k(l-A(k)) 



(m + B(k)) 



1/2 ' 



(8.1) 



pA(k) 



+ 



(k 2 (l- A(k)) 2 + (m + B(k)) 2 ) 



1/2 ' 



(8.2) 
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where we have defined the potential V as 



V(p) = — 8iia^-, p = \p\, V(x) = ar, r — |x|. (8.3) 



P 4 ' 

In the chiral limit m = they become: 

2J ( 2?r ) (fc 2 (l-A(A0) 2 + 5(fc) 2 ) 7 

1 f d 3 k , . k(l - , s 

pAp = + -/ — — Fp-k- ^ ^ 8.5 

These equations were extensively studied in the literature; they result from the as- 
sumption that confinement is due to the exchange of Coulomb-gluons [3U] I3T] 132] 135] 133] . 
Using variational methods, it has been shown that the perturbative vacuum state is un- 
stable under quark-antiquark pair creation. The existence of a new stable vacuum state 
is ensured by the existence of a nontrivial solution to Eqs. (|8.4|H|8.5|) . We summarize 
below the main results that have been obtained and outline some salient features related 
to physical aspects of the problem. 

Equations (|8.4|) - (|8.5|) are solved by decomposing the functions B and p(l — A) along 
polar combinations, by introducing an angle and a modulus: 

B{p) = E(p) sin<£>(jo), p(l — A(p)) — E(p) cos <p(p). (8.6) 

Equations (|8.4j) - (|8.5[) can be rewritten in the form: 

1 r d 3 k 

E(p)smp(p) = --jj^V(p-k)smp(k), (8.7) 

1 r d 3 k - p-k 

E(p)cosip(p) = p- - J -^-^V(p- k)p.k cos ip(k), p = -, k = -. (8.8) 

The latter equations in turn can be recombined to decouple the function tp from E: 

1 r d 3 k 

psmip(p) = — / - — r^y(p — k) sm<p(k) cos ip(p) — p.kcos ip(k) siny?(p) , (8.9) 

2 J [2ix) A L - 1 

1 r d 3 k 

E{p) = pcos(f(p) — -J j^-^V(p — k) smip(k) sin (p(p) + p. k cos (p(k) cos (p(p) 

(8.10) 

From Eq. (J8.9|) one deduces that the function if is an infrared finite quantity. From Eq. 
(|8.10|) . after expanding in the integrand k-dependent terms around p, one deduces that E 
is infrared singular and its singularity is represented by the integral —\ J d 3 kV(k)/(2n) 3 , 
a property mentioned and utilized at the end of Sec. U\ for the checking of the infrared 
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finiteness of the wave equation (J7.12J1 . (We use, however, dimensional regularization for 
analytic calculations throughout this work.) 

The function sin </?(/?) wil be identified later with the Goldstone boson wave function. 
Therefore, it should be a normalizable and presumably nodeless function. An analysis 
of the above equations shows that the function ip behaves at infinity as p~ 5 and tends 
at the origin to ir/2. The solution that is found is indeed a monotonically decreasing 
function with the above properties. On the other hand, due to the infrared singularity of 
the potential V, the energy function E vanishes at some finite value po of p and becomes 
negative for p < po; it tends to a finite negative value at the origin and behaves as p at 
infinity. The functions B and p(l — A) vanish simultaneously at po; they are negative for 
p < po, with the limiting values B(0) = E(0) and lim p=0 p(l — A) = 0; asymptotically, B 
behaves as p~ A and A as p~ 2 . 

An order parameter for chiral symmetry breaking is the quark condensate < ipip >. 
For a given type of quark, say it, it can be defined as minus the trace (in color and Dirac 
spinor spaces) of the quark propagator at the origin in x-space: 

< uu >= -tT csp S(x) = -N c ( -^-rtr sp S(p). (8.11) 

Using definition ()7.6|) and decomposition ()7.8|) . one finds: 

< uu>=-2N c J -0^ sin <p(p). (8.12) 

From Eq. (|8.7|) one also sees that the quark condensate enters in the asymptotic behavior 
of the mass term of the quark propagator in the chiral limit: 

psmip(p) ~ -— — . (8.13) 

This relation is the analog of that obtained in perturbation theory from operator product 
expansion and the renormalization group analysis 

We next turn to the bound state equation (|7.12|) . A consistent study of the chiral limit 
should be done starting from the situation where the quark masses are different from zero 
and then taking the limits mi = m 2 = 0. In the present case, we directly consider the 
equal mass case ni\ = ni2 = m, corresponding to the light quark sectors u and d without 
isospin breaking. 

In general, the resolution of Eq. ()7.12j) proceeds by first decomposing the wave function 
ij) along a basis of 2 x 2 matrices spanned by the matrices 7^ and 75. In the following, we 
consider the rest frame of the bound state (P = 0). One has the decomposition 

i> = ^1 + 70^2 + 75^3 + 7075^4, (8.14) 
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where the functions tjj a , a = 1, ... ,4 are themeselves 2x2 matrices on which act the 
spin Pauli matrices er. The quantum numbers of the state are usually defined with 
respect to the components that survive in the nonrelativistic limit; these are ^3 and i/j^, 
which have the same quantum numbers. For the present problem, we are considering the 
sector characterized by the following quantum numbers: total spin s — 0, orbital angular 
momentum i = 0, total angular momentum j = 0. We adopt the rule of writing the 
energy sign operators present in the potentials [Eqs. ()6.20|) and (j7.10|) - (j7.11|) ] as well as 
the factors ((1 — A)/E)p t on the utmost left. The doubly transverse momentum operator 
p* [Eq. (J6.17)) ] annihilates S"-states; it commutes with x 2 and with all scalar operators 
that act on S"-states in x-space. In the equal-mass case, the transverse potentials and 
A l 2 are opposite to each other; defining 



A< = ^^P*A = -A*, (8.15) 

using definition ()8.Hjl . introducing the spin operators S\ and S2 of particles 1 and 2, and 
removing indices 1 and 2 from various factors which are equal, we obtain the following 
four equations: 

Po^i + 2(1 - A)( Sl - s 2 ).p(l + ^V)ij 3 = 0, (8.16) 
P01P2 = 0, (8.17) 
P ^3 - 2(m + B)(l+ ^V)^ + 2(1 - A)(si - s 2 ).p(l + ^V)^ 

-2 (1 ~ A) (si - s 2 ).p^i = 0, (8.18) 
hi 

P ^4 - 2(m + B)(l+ ^V)^ = 0. (8.19) 

We have neglected in Eq. ()8.18)1 the internal L 2 dependent parts, present in the potentials 
(Ail + A 2 l) and A* acting on ipi, which is a P-state, and used the limits (j6.2Hj) - (Jfi.24j) : 
these neglected parts do not seem to play a crucial role in the following calculations; with 
the decomposition ()8.15|) . A = ar/6. 

From Eqs ()8.16|) and ()8.19|) one deduces the relation: 

(m + P)^i + (1 - A)(S! - s 2 ).p^ 4 = 0. (8.20) 

This equation, which is valid in general for m 7^ and Pq 7^ 0, should also be satisfied 
in the limit m — > and for Pq — > 0, if the Goldstone boson properties are expected to 
vary smoothly under explicit chiral symmetry breaking. Then, the ground state solution 
of the above equations in the limit m = with zero energy is: 

4>i = 4>2 = ^4 = 0, ^3 7^ 0, (8.21) 
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with ifj 3 satisfying the equation 

1 r d 3 k 

E( P )Up) = — 2 j ^s^(P - k )^)- (8.22) 

This equation is similar in form to Eq. (|8.7|) . indicating that ip^ is proportional to the 
self-energy function sin ip: 

Va(p) =Csin^(p), (8.23) 

with C a constant. 

In order to calculate the pion decay constant F n , one must again consider the case 
where m^fl. The pion decay constant is defined as: 

< 0|d(0) 7M 75^(0)|7T + (P) >= iV2F^. (8.24) 

This is related to the component ip^. 

k ;Mk) = ^F n P . (8.25) 



(2tt) 3 ^ v ' 4 

The normalization condition ()6.25|) yields, after elimination of ip\ through Eq. ()8.20|) and 
replacement of ^3 by its expression ()8.23|) : 

f d 3 k 

8C J j^Mk) = 2P N C . (8.26) 

Using Eq. (|8.25|) . one finds: 

CF n = (8.27) 

This result is actually a consequence of the Ward-Takahashi identity relative to the Green 
function < Tu(x)'-f fl '~f5d(x)id(0)'~f5u(0) >, which implies (in the chiral limit) the equation 
\/2F w < 0|i<i(0)75u(0)|7r + >= 2 < uu >; this is equivalent to the relation 

d 3 k , . <uu> . . 

Mk) = -—>=—. (8.28) 



(2vr) 3r " v ' 2y/2F„' 

Using Eq. (jHHJ), one obtains Eq. 1|P7) . 

The normalization condition (|8.26|) can also be analyzed by using Eq. ()8.19j) . Inte- 
grating that equation with respect to p and expanding in the second term all quantities 
with respect to m to first-order in m and then using the integral equations satisfied by 
the first-order self-energy functions, one ends up with the relation 

d 3 k , 2m f d 3 k , 

i1>4 = -5- 77^3, (8-29) 



(2vr) 3r " P J (2tt 



35 



which is equivalent to the well-known relation of Gell-Mann, Oakes and Renner 
m^F% = —2m <uu>. 

In order to be able to calculate F n , one needs to know the function This information 
comes from Eq. (j8.18|) . Eliminating ipi through Eq. ()8.2()j) and noticing that ip4 iself is 
proprtional to Po, one can take in all factors multiplying ip4 the chiral limit. Defining 

^ 4 = P ^ 4 = P ^3, (8.30) 

the function g(p) then satisfies the integral equation 

E(p)g(p) = sin (p(p) 

1 { d 3 k r * 1 

- 2 J v (p - k ) [ sin sin <p(p) + p- k cos cos <p(p)\ y(k) 

1 r d 3 k 

+ -cos^(p) J — ^T/(p-k)(p-k).kcos^(%(fc). (8.31) 

Except for the last term, which comes from the contribution of the rotational motion of the 
flux tube, this equation is the same as a corresponding one derived in Refs. |40| I42j. The 
infrared singular part of E cancels the singular part of the first integral and the function 
g appears as infrared finite. This is also seen more explicitly by using Eqs. ()8.9j) - (j8.1()j) : 
after a few algebraic manipulations one can cast Eq. ()8.3H) into the form 



pg{p) cos(/?(p) = sm. ip{p) 
1 f d 3 k 

+ 2 J J^3 V ( p ~ k ) { sini P( k ) sin<^(p) + p.kcos^(fc) cosp(p)] (g(p) - g(k)) 
1 r d 3 k 

+ -cos^) J y(p-k)(p-k).kco S ^(%(A;). (8.32) 

which is free of infrared singularities. The function g has the same asymptotic behavior 
as sin ip and turns out to be nodeless. Integrating this equation with respect to p yields 
the relation 

; kg(k) cos (f(k) = [ sin (p(k) 



(2tt) 3 y j rv ' J (2tt) 

r d^p 1 r d^ k * 

+ J (^)3^ cos ^) J J2k)Z V ^ P ~ k ^ P - k )- kcos ^(%( fc )> ( 8 - 33 ) 

which can be considered as a consistency check for various approximations and numerical 
calculations. The resolution of Eqs. (j8.31jl or ()8.32jl provides the function g and hence 

A rough evaluation of the orders of magnitude of the quark condensate and of the 
pion decay constant can be done using analytic approximations. Taking into account the 
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properties of the function ip at infinity and near the origin and its nodeless character, the 
following approximations can be used inside the integrals: 

• (b\ 1 ,,s k f 1 y/E-1 \ 

Smip{k) ~ (l + V/Qxr))*/*' C ° SIP{ ' ~ ^[(TTfcTWF + (l + ky(ba))^r 

(8.34) 

where b is a constant to be determined. For the numerical applications we use a = 0.18 
GeV 2 , \fa = 424 MeV. With these approximations, the integrals in Eqs. (|8.7j ) -([8.8 |) 
can be evaluated analytically (with dimensional regularization) yielding hypergeometric 
functions. The parameter b is then determined by the requirement that the left-hand 
sides have the same zero. We find b = 0.31. The evaluation of the integral in Eq. (J8.12)) 
(with approximation (J8.34J) and N c = 3) gives < uu >= —(115 MeV) 3 , in agreement with 
the results found in Refs. USUI, < uu >~ -(100 MeV) 3 . 

Since the function g has similar asymptotic properties as simp, we approximate it as: 

g(p) - -7j=siny>(p), (8.35) 

where 7 is a constant; it is determined from Eq. ()8.33|) . At a first stage we neglect the last 
integral; we find 7 = 0.52. Equations 1ET2EI) . (1^271) and ((ODD then S ive K = 16 MeV, 
with the same order of magnitude as the results found in Refs. |l2] El] . F n ~ 11 MeV. 
Consideration of the last integral in Eq. (J8.33j) increases the value of 7 by an amount of 
16% and brings F n to 18 MeV; it does not therefore change its order of magnitude. The 
experimental value of F n is nearly 93 MeV. The quark condensate < uu > is not a directly 
measurable quantity; QCD sum rules [17] give the prediction < mm >~ —(225 MeV) 3 at 
the scale of 1 GeV. 

In Ref. jl^ , numerical calculations have been done with the harmonic oscillator 
potential, which makes it difficult to compare predictions. In Refs. |4*Ul I43j . rather large 
values of F n and < uu > are presented, but these are obtained by fitting parameters in 
order to fix the value of the constituent quark mass at 300 MeV, corresponding to big 
values of the string tension (or of its equivalents). 

The relative smallness (about one order of magnitude) of the theoretical predictions 
for <uu> and F v found here and in Refs. O El] seems to indicate that short-distance 
forces should have sizable contributions in these quantities; this is corroborated by the 
facts that, first, they act in momentum space with a potential having the same sign as the 
confining one, and second, they induce much slower asymptotic behavior to the function 
tp. A definite conclusion about the quantitative aspects of chiral symmetry breaking could 
be reached only when both kinds of forces, long-range and short-range, are considered. 
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9 Properties of the bound state spectrum 



We turn in this section to a study of the main qualitative properties of the bound state 
spectrum that emerge from the wave equation, leaving to a separate work a more quan- 
titative study of it, in particular when short- distance effects are incorporated. 



9.1 One-particle limit 

We begin with the case when one of the particles becomes infinitely massive; this is 

achieved by taking one of the masses, m 2 , say, to infinity. Defining P = m 2 + pio, 
x = X! — x 2 , p = (px — p 2 )/2 = pi, the wave equation ()7.12|) reduces to a Dirac type 
equation for the quark in the presence of a vector static potential: 

p w - h - A + 7io7-A]^ = 0, (9.1) 

where the potentials, in the classical limit, are: 



A = ff H\/^arcsiny^-(l- e ( Pl0 ))^(l-^l--^)^, (9.2) 



ar f p*x E 2 x 2 f £ 2 x 2 /L 2 / L 2 1 

= iy-^Wi^^fei l -E^y (9 - 3) 

[r = |x|, p* is defined in Eq. f!6.17|) . L is the orbital angular momentum operator, hi, 
Ei and e(pio) are defined in Eqs. (j7.9|) - (j7.11|) .] They represent the energy-momentum of 
a straight segment of length r with linear energy density a turning around the point x 2 
(the position of the heavy antiquark). 

Taking further the limit of a heavy quark, the potentials become to order 1/c 2 : 

Ao = *r(l + -^), A = ^. (9.4) 
v bmfx. z/ 3 mi 

9.2 Nonrelativistic limit 

Next, we consider the case of two heavy quarks. The potential Ai^ [Eq. (Jfi.2()|) ] becomes 
(in the cm. frame) to order 1/c 2 : 



r l 1/ mi w 1 1 1 \L 2 

A w = ar 



■)(- 



2 ' 6^mi + m 2 ^m 2 m\ m^'x 2 



1 / mim 2 \2 / 1 1 \ L 



2_ 



m 1 +m 2 / v mf m 2 'x 2 
1 

3mi 3m 2 



(9.5) 
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the definitions of x and p being the same as in Eqs. (J6.3|) . A2 is obtained from A± by 
exchange of indices 1 and 2 and replacement of p by —p. The self-energy £ [Eqs. ([7.8)1 . 
()8.1)l - ([8.2j) ] behaves for large m as: 

£(p,m) = -— + 0(l/m 2 ). (9.7) 

After having subtracted the masses and made a few changes of function to reach an 
explicitly hermitian form, the nonrelativistic hamiltonian, to order 1/c 2 takes the following 
form (in the cm. frame): 

# = 7r + ar — + — )- o(— 3 + — s)(p ) +^r(— 2 + — 2)- 

2\i tx x mi m 2 y 8 v mf my 4 v mf my r 

a / 1 1 1 \ /T 9 cr /L.si L.s 2 

-2r(— + — )( L + 2h + ;rHr + 

2(7/ 1 1 x 2(7/1 1 \ 

(—5 )L.Si (—5 )L.s 2 . 9.8 

3r v mf 2m\m2 J 3r^m^ 2m\m2' 

[p, = mim 2 /(mi+m2), Si and s 2 are the spin operators of the quark and of the antiquark.] 
Several remarks can be made at this stage. First, the hamiltonian is independent of spin- 
spin interactions; this is already evident from the wave equation (|7.12|) where no direct 
interactions involving the spacelike 7-matrices of both quarks exist. The absence of long- 
range spin-spin interactions is compatible with experimental data. Second, we notice the 
presence of purely orbital angular momentum dependent pieces (proportional to L 2 ), the 
origin of which is related to the contribution to the rotational motion of the system of 
the moment of inertia of the flux tube, represented by the straight segment joining the 
quark to the antiquark. The corresponding centrifugal energy produces a global plus sign 
in front of those terms; the minus sign results from the additional contributions of the 
momentum of the flux tube, which also couples, through the spacelike potentials Ai and 
A 2 to the quarks [Eq. ()9.6|) ]. Those terms were also obtained in Refs. [TH] and jMJES]- 
Third, we have two kinds of spin-orbit term. The first, appearing after the term in L 2 , 
comes from the contribution of a conventional timelike vector interaction represented by 
the potential or, which is the dominant part of the combination A\q + A20 [Eq. ([9.5)1 ] . 
The second type of contribution, provided by the last two terms of Eq. ([9.8)1 . comes from 
the contributions of the direct interactions of the momentum of the flux tube with the 
quarks, represented by the spacelike potentials Ai and A 2 [Eq. ([9.6)1 ] . The latter terms 
induce negative signs to the spin-orbit couplings, in opposite direction to the former one, a 
feature which is also observed on phenomenological grounds for the large-distance effects 
in fine splitting. 

Adopting the notations of Ref. ^Hj for the potentials corresponding to spin-orbit (Vi 
and V2) and to spin-spin interactions (V3 and V4), we have: 

7i = ~<rr, V 2 = +^ar, V 3 = 0, V 4 = 0. (9.9) 
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Designating by V the nonrelativistic confining potential [Eq. ([8.3)1 ] . the potentials V, V\ 
and Vi satisfy the Gromes relation 

V + V 1 -V 2 = 0, (9.10) 

which is a consequence of Lorentz covariance. 

The expressions of potentials V\ and V 2 do not satisfy, however, Buchmiiller's con- 
jecture about the spin-orbit terms [2H]> according to which confinement is due to a pure 
color electric field in the co-moving frame of the two quarks, thus reducing the spin-orbit 
potential to a Thomas precession term and entailing V\ = —or and V 2 = 0. Leaving 
aside the question of a phenomenological determination of these potentials, we notice 
that lattice calculations [4*0] 150] . which favor the negative sign of Vi, do not, however, 
clearly distinguish between the two possibilities above, due to the existing uncertainties. 
In order to clarify the structure of the spin-orbit potentials, it is necessary to investigate 
the contributions of higher-order terms not taken into account in the potentials (|fi.20|) (of 
the type of those considered for the extraction of the self-energy terms). Another inde- 
pendent check consists in calculating the field correlators appearing in the nonrelativistic 
expansion of Ref . [IB] . 

Reviews about bound state problems of quarks can be found in Refs. [HI] and [5*2"] . 

9.3 Regge trajectories 

Finally, we consider the high energy behavior of the mass spectrum for ultrarelativistic 
systems. It has been known for a long time that the linear confining potential of the 
static case produces in the ultrarelativistic limit of massless quarks, through the Salpeter 
equation, linear Regge trajectories. The inclusion of a flux tube, represented by a straight 
string, modifies the relationship of the angular momentum with the total mass of the 
system and increases the slope a' of the Regge trajectories by an amount of 15-20%, 
enforcing the classical relation a' = (2ira)^ 1 with the string tension a [o"3*] I35j. 

We have checked these properties on Eq. ([7.12J1 . by solving it in an approximate 
way that preserves its main qualitative features. The approximation that we use is the 
Breit approximation j20|, which consists in transforming the wave equation into a local 
equation in x-space. The main lines of the approximation are presented in appendix IU1 
The resulting Regge trajectories for the tt and p families are presented in Fig. El where 
for comparison we have also presented the trajectories of the linear confining potential 
corresponding to the situation where A 10 + A 2 o = or and A x = A 2 = in Eq. ([7. 12)) . 
One verifies, first, the linearity of the trajectories and, second, the increase of the slopes 
when the flux tube is present. Similar trajectories and behaviors are also obtained with 
the cio and a\ families. 
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There are uncertainties of a few percent in our results coming from the way of cal- 
culating the mean values of operators present in the flux tube functions, depending on 
whether one evaluates them on the full function ijj or only on the basic component ^3 
[Eq. (j8.14j) ] and whether one uses for the latter a free norm or the norm relative to the 
Pauli-Schrodinger equation that it satisfies. 
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Figure 6: Regge trajectories of the it and p families and their first daughter trajectories. 
Full lines correspond to the flux tube potential, dotted lines to the linear potential. Inputs: 
m x = m2 = 0, a = 0.18 GeV 2 . 

On phenomenological grounds, one actually determines the value of the string tension 
a from the experimental Regge trajectories using the string theory relation a' = (27rcr) _1 . 
From the slope of the 7r-family trajectory one obtains a = 0.22 GeV 2 , while from the p- 
family trajectory one obtains a = 0.18 GeV 2 A precise comparison of our results with 
those data necessitates a more elaborate resolution of the wave equation and presumably 
the inclusion of the short-distance potentials. 
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10 Summary and concluding remarks 



We have investigated, in the large N c limit, the large-distance dynamics of QCD within 
the quarkonium system through the saturation of Wilson loop averages by minimal sur- 
faces. The dynamics is described by a wave equation of the Breit-Salpeter type, in which 
the interaction potentials are provided by the energy-momentum vector of the straight 
segment joining the quark to the antiquark and carrying a constant linear energy density, 
equal to the string tension; the latter represents the effective contribution of the color flux 
tube of the quarkonium system. In the static case, confinement is realized by the usual 
linear potential, while in the general case of moving quarks additional contributions come 
from the moments of inertia of the straight segment . 

Taking into account the self-energy parts of the quark propagators, it was shown 
that chiral symmetry is spontaneously broken with a mechanism identical to that of the 
exchange of one Coulomb-gluon. 

In the nonrelativistic limit, long range spin-spin potentials are absent, while the mo- 
ments of inertia of the flux tube bring contributions to the orbital angular momentum 
and spin-orbit potentials with negative signs, in opposite direction, for the latter case, to 
the spin-orbit term of the pure timelike vector potential. 

In the ultrarelativistic limit, the mass spectrum displays linear Regge trajectories, the 
slopes of which tend to satisfy the classical relation with the string tension obtained in 
straight string theories. 

The potentials that have been isolated in the present work do not represent the com- 
plete interaction kernel of the wave equation, but only the large- distance dominant con- 
tributions to it. The determination of the higher-order terms necessitates further study. 
Also, mathematical assumptions made for neglecting certain terms or approximating oth- 
ers in the large separation time limit need to be analyzed in more detail. 

The formalism that has been developed here, centering the Wilson loop representa- 
tion on minimal surfaces, ignores the contributions of short-distance effects provided by 
perturbation theory. A complete resolution of QCD requires the presence of both large- 
and short- distance effects. A corresponding resolution of the loop equations is far from 
being trivial, although approximate solutions might still provide a practical framework 
for a quantitative investigation of the spectroscopic properties of the theory. 
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A Properties of minimal surfaces 



An important quantity in loop dynamics is the second-order functional derivative of the 
area at two different points on the contour. For definiteness we consider the two points on 
the line r = (line X1X2 of the contour C of Fig. For a minimal surface, the variations 
of the contour maintain at all stages the minimality property of the area; therefore the 
second term of the right-hand side of Eq. ()3.7|) and all its variations are null and can be 
ignored. The second-order functional derivative of the area is then: 

5 SA 8 r 

(S^p - S^a)— 8(a - a') / da^(y)F(y - x(a, 0)) 



5xP(a>,0) 5x"(<r,0) y ^ a vp ^p~™> da 
+ (W>> °) - <Wfa 0)) (^x'V, 0) - 5„ p x'»(a', 0))f(x(<t, 0) - x(a', 0) 



+ (5, a x'»(a,0)-5„ a x'»(a,0))5(a-a') [ da^(y)F(y - x(a, 0)) 



-(5 lxa x' u (a,0)-8 l/a x^(a,0)) J da" dr 



dxP(a,0) 
,_„ 5y\*",T") 
5xP(a',0) 



(y'if - + (?/V - z/VKa + (yV - ^f{v - o)). 

(A.l) 

The evaluation of the the last integral can be done by expanding 5y x /5x /3 (a' , 0) about 
the point x(a, 0) that appears in F: 



X(„n 

+ ..., (A.2) 

r=0 



5y x (a",r") 6y x (a,0) , d <V(a,0) „ d 5y x (a,r 

(a — a) — „ „, , -t + r 



5xP(a',0) 5xP(a',0) v ' da 5xP(a', 0) dr5xP(a',0) 

where the dots stand for terms that do not contribute in the limit a = 0. One notices that 
the first term in the right-hand side of Eq. ()A.2jl is a derivative that lies on the contour 
and hence is equal to 5 ^5 (a — a'): 

5y x (a,0) _ Sx x (a,0) 



5 px 5{a-a'). (A.3) 



It is also convenient, to complete the calculations, to stick to the following diagonal 
constant metric: 



t x /2 



x 2 



x'.x = 0, x 2 = const., x 2 = const., A = \j (A.4) 

The result is: 

2 ±S(a-a f ) I da af3 (y)F(y - x(a, 0)) 



5xf 3 (a',0)5x a (a,0) da 

+2(5 a px'(a, 0).x'(a\ 0) - x ,f> (a, 0)x /a (a\ 0))f(x(ct, 0) - x(a\ 0)) 
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'x' 2 d Sy x (a,r) 



x 2 dr 6xP(<7', 0) 



T = 



x' x (a,0)x' a (a 1 0) x A (a,0)x>,0) A 

x*(a,0) x2(a,0) > ^ 



The two points x(<r, 0) and x(a', 0) lying on the boundary, the function F\x(a : 0)— x(a' : 0) 
is equivalent, in the limit a = 0, to 2a% / 3 ./ 2 ( H er — O an d thus cancels the singular term 
of the third line. (We assume here that the line under consideration does not have 
self-intersections neither displays backtracking.) There remains to evaluate the quantity 
Jr ifohy"ol ' This i n turn requires a more detailed knowledge of the derivative |^ . 
That quantity does not have a compact expression for general contours. Studies of its 
main properties can be found in Refs. jSJ and In the formulation that we present 
below, we follow partly the covariant approach adopted in Ref. j2H] and put emphasis on 
the aspects that will be needed in the present work. 

Let us first point out the fact that a variation Sx on the contour that lies in the tangent 
plane to the minimal surface does not modify the internal part of the minimal surface but 
only adds to it a new small piece on its boundary within the tangent plane; the variation 
thus amounts to a simple extension of the boundary on the same surface. In that case Sy 
should lie in the tangent plane to the minimal surface at y, being a linear combination of 
y' and y and corresponding to a reparametrization of the surface. In this respect, it can 
be checked that if 5y x contains terms proportional to y' x or to y x then the contributions 
of these in Eq. (|A.1J) are identically zero. Therefore, the nontrivial contributions to Eq. 
(|A.1J) come only from the transverse variations of Sy with respect to the surface, which 
are themselves due to transverse variations Sx on the contour with respect to that surface. 
We shall therefore concentrate on these kinds of variation. 

To proceed further and to exhibit some general features, we shall use for a while 
covariant notations. The parameters of the surface will be denoted £ a , a = 1,2, and the 
metric g ab : 

9ab=g^g^, a, b = 1,2. (A.6) 

The entries of the inverse of g are denoted g ab and raising and lowering of indices are done 
with these tensors. The determinant of g is denoted \g\. We shall not meet below singular 
quantities and therefore, as long as not needed, we shall not use regulated expressions. 
The area is then A = J d 2 ^\g\ 1 / 2 and the equation of the minimal surface is: 

1/2 d a (\g\ 1/2 g ab d b y l *(O)=0. (A.7) 



\g\ 

The determinant \g\ remains constant on the minimal surface. 

Next, we introduce variations Sy on the minimal surface, including the boundaries. 
The equation of the new minimal surface is obtained by replacing in Eq. (|A.7|) y by y + Sy. 
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The linearized equation in 5y is: 

A, u 8f = \g\-^d a (k| 1/2 M;W 

where one has defined 



M, 



5»u9 ab + try,d b y u - d a y u d b y, - g ab d c y v d c y,. 



ab ■ 



fit' 



The operator A^ v is transverse: it annihilates any tangential variation Sy L 

A. v u c d c y» = 0. 



(A.8) 

(A.9) 
uj c d c y u : 

(A.10) 



One next introduces the Green function of the operator A^ subject to the boundary 
condition that it vanishes along the contour C: 



A^G^iCv) = \gr 1/2 (s, P -Q, P (Oy^-v) 



(All) 



where is the projector on the tangent plane of the initial minimal surface at ?/(£): 

Q,M) = d a y^)g ab d b y P (0- (A.12) 

Using then the Green identity, the symmetry property Mj%, = M^" the Stokes theorem 
and the transversality of 5y, one obtains for the transverse variations 5y the following 
successive expressions: 



'^ a [^(6l3l 1/2 M p t(a^(^) - H 1/2 mJ^i/,(0)g w (^)' 

/ d^ca\g\ 1/2 g ab 5yv(Od b G ull (^ v ) 

J C 

dXn a (X) 5y v (£ (A) ) d a G V)Ji (£ ( A) , 77) , (A.13) 



c 



where the contour C has been parametrized with the parameter A and where n is an 
orthogonal vector to the contour in the tangent plane to the minimal surface, 



n a (X) = g ab \g\ 1/2 e b( 



The functional derivative of y^f]) with respect to x u (l;(X')) then becomes: 

8yM 



M£(A')) 



n a (X')d a G^, V ) 



S=«(A') 



(A.14) 



(A.15) 
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Finally, the derivative of along n is: 
%(»7) 



n 



'<M£(A')) 



A consequence of Eq. (jA.16|) is the symmetry property of its right-hand side with 
respect to the variables £ and t] for A ^ A', a feature that is not evident in the left-hand 
side. That property is intimately related to the symmetry property of the second-order 
variation of the minimal surface with respect to the order of operations. To exhibit this 
feature more explicitly we return to the (a, r) parametrization that we were using above. 
The derivation n a d a is equivalent in the limit r = to the derivation — v/^ J~- Equation 
(|A.15|) becomes, using the metric (jA.4j) : 



Sy x (a, t) 



5xP(a',0)) V x 2 dr' 
which in turn yields for Eq. (jA.5|) the expression 
5 5 A d 



7 G px (a',T';a,r) (A.17) 



5xP(a',0) 5x a (a,0) 

J2 



2—5(a-a') J da^{y)F{y - x(a, 0)) 



x d d 



S\n — 



x' x (a,0)x' a (a,0) x x {a,0)x a (a,0) 



r'=0,r=0\ X' 2 X 2 

(A.18) 



The symmetry property of the second-order derivative of the minimal area can now be 
discussed for the case a ^ o' . The Green function can be decomposed along symmetric 
Lorentz tensors each multiplied with a scalar function. Because of the boundary conditions 
(jA.llj) the two derivatives and must simultaneously act on at least one of the scalar 
functions in order not to yield zero. In that case the product of the transverse projecter 
present in the last term of Eq. ()A.18|) with the tensor factors of G$\ gives a symmetric 
result, which ensures the commutativity of the two functional derivatives when acting 
at two different points of the contour. That symmetry property is evident when the 
functional derivatives act on the Wilson loop average; then the two operators manifestly 
commute. This study will be completed below by also including the case of coinciding 
points. 

Concrete calculations can be done when the Green function is specified more explicitly. 
To this aim, one can try to construct the Green function by means of an iterative procedure 
with respect to powers and derivatives of the curvature of the minimal surface: 



^(£,7/) = ££«(£,/,), (A.19) 



i=0 

where the index % is related to the power or the order of the derivative of the curvature. 
The lowest-order term of this iteration would correspond to the case of a plane with zero 
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curvature. The solution to that problem can be explicitly constructed. Using the metric 
7m. one can use the ansatz 



G{ p\Wy\v, r ) = Up\(z'.y'z.y - z'.yz.y') 

-z-y{z' x y' p - z' p y' x ) + z'.y(z x y'p - z p y' x ) 

-z'.y'{zxyp - zpyx) + z.y'(z' x y )3 - z^ A )|-^/i {0) (a', r; cr, r), 

(A.20) 

where we have defined z = x(a',T f ), y = x(a,r) and kept in the decomposition of Gp X 
tangential components in order to satisfy at the end the boundary condition ()A.3j) . The 
function is assumed to satisfy the inhomogeneous Laplace equation 

A {0 h(°\Z; V ) = A iv) hW(Z; V ) = 5 2 (Z- V ), £=(a',r'), ij = (ff,r), (A.21) 
1 d 2 d 2 

A « = w + V- (A - 22) 

and the boundary conditions 

h(^v)\^ c = h^v)l £C = 0. (A.23) 

The coordinates z and y represent points on the background minimal surface and with 
parametrization (jA.4|) satisfy the Laplace equations 

A (0 ^(0 = 0, \ri)yM = 0- (A.24) 

Inserting G$ x in Eqs. (|A.11|) one can reconstruct G^ x by appropriately decomposing it 
along Lorentz tensors as in Eq. (jA.20|) and determining its Lorentz scalar parts through 
inhomogeneous Laplace equations involving as inhomogeneous terms and its first- 
order derivatives. It should be noticed that the inhomogeneous part of Eqs. (jA.ll|) is 
already saturated by the function G^ x (the factor | 1 1 / 2 being absorbed in the definition 
of the Laplace operator (| A. 22(1 ) and therefore the higher-order functions do no longer 
have contributions to the delta-functions. The Lorentz tensor parts of the function G^ x 
contain at least one second-order derivative of y or z and therefore are of first-order in the 
curvature of the surface. The procedure can be continued to higher orders; the function 

(2) 

Gp X will be quadratic in the second-order derivatives of y and z, or will contain third-order 
derivatives, and so forth. On technical grounds, it is advantageous to use also a tensor 
notation for the parameters of the surface and to simultaneously decompose the higher- 
order functions along tensors (defined essentially by the derivative operators) in parameter 
space. We shall not develop here the above procedure, for we shall not need in the present 
paper the explicit expressions of the functions G*\ (i > 1). In many instances, the sole 
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knowledge of the function is sufficient to obtain or to check the main properties of 
the system under study. This is in particular the case where the background minimal 
surface from which deformations are calculated is itself a plane. 

Let us now concentrate on the function Gffj. We consider a four-sided contour C as 
in Fig. with the additional restriction that the four sides correspond to the lines of 
equations r = 0, a — 0, r = 1, a — 1. The explicit expression of the function h^°\ 
satisfying the boundary conditions (jA.23|) is then: 

,(n\,,, s 4 ^ sinfnvro"') sin(n7ro") sin(m7rr') sin(m7rr) 

h^(a',r';a,r) = £ — ^ >- ) ' \ ^ >-. A.25 

7T 2 n ^ =1 (m 2 A + n 2 /X) 

Calculating, in view of Eq. (|A.17J) . the derivative of with respect to r' at r' = 0, one 
obtains a series in m that can be summed into a hyperbolic function: 



,%( ) 00 smimra" ) smiwrc) smn I wri l — t)/a ) 
/(°)(a',0;a,r)^-A^ = 2£ — -J n \ 11 • (A.26) 



dr' 



sm(mra') sin(n7rcr) sinh (rnr(l 



t'=o ^ sinh(wr/A) 



One checks that satisfies, with and through Eq. ([A.17J) . the boundary condition 

dH). 

We next evaluate the derivative of P-p with respect to r at r = 0. The derivative 
acts on f( QS> and on the Lorentz tensor part of The latter action yields in Eq. (|A.18|) 
the finite part of x' 1 ' J da au (y)F(y — x(a, 0)) multiplied by 5(a — a') coming from the 
restriction of on the contour. One finds for the second-order functional derivative of 
the area the result: 

5 5 A 3 



2— 5(a - a') [ da aP (u)F(u - y) 
da J 



+25(a - o')y' v ^- n I da au (u)F(u - y) 



5xP(a',0)5x a (a,0) 

dy 13 J f p 
5p a (z'.y'z.y - z'.yz.y') 

1 0/<°V.O;ff,r) 



y,2 




y 2 



{z-yZaVp + z '-y'z a yp - z'.yz a y' p - z.y'z'Jp) 



y' 2 y 2 dr 



T=0 

(A.27) 



where y = x(a, 0), z = x(a',0) and the subscript "fp" means finite part. The last term 
of Eq. (|A.27|) is invariant under the change of the order of derivations, as can be seen 
from the tensor terms and with the use of the defining equations (fA~25jl - (fA~26|) of / (0) , 
reflecting the general property emphasized after Eq. (jA.18j) . The symmetry property of 
the contact terms will be established shortly. 
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The singular part of 9 ^ ) ^ r '°'°''' r ^ 
from Eq. (lA~26j) . It is: 

dfM(a>,0;a,T) 



T = 



when (a — a') approaches zero can be calculated 



Or 



1 / x 



=0 (a-a')-*0 7T V X 




a — a 



A-2 



0((a-a')°). 



(A.28) 



It yields for the contracted second-order derivative 
terms do not contribute, the behavior 

S 2 A 



S 2 A 



5x a (a' ,0)Sx°> {a,0) ' 



in which the contact 



a - a')~ + 0((a - a 



/\0\ 



(A.29) 



5x a (a',0)dx a (a,0) (a-a')^o n 

It is in agreement with the leading part of the short- distance behavior obtained in Ref. 
|26j . The contracted second-order derivative contains also a logarithmic singularity 
which is obtained in the present approach by considering the contributions of the next 
two nonleading terms in the curvature, i.e., the functions Gfy, i — 1,2 [Eq. (|A.19|) ]. 

We now check the symmetry property between two successive functional derivatives. 
Let us first consider the problem in the Wilson loop case. The commutator of two func- 
tional derivatives on the contour is, using Eqs. (|2.3|) and (J2.1UD : 

S - i *—\ W(C) = 

d d 
ig— 5(a - a'){F Pa {x{o))) w - ig—5(a - a'){F a p(x(a'))) w 

+ig5(a - a')x ,v (a)((ypF ua (x(a)) - V a F u/3 (x(a)))) w . (A.30) 

Expanding in F a p(x(a')) a' about a and using in the last expression the Bianchi identity 
one finds: 



5 



5 



5 



8xP(a')8x a (a) ox a (a)8xP(a') 



W{C) 







^da 
d_ 

da 







^(S(cr'-a)(a'-a))^(F Pa (x(a))) 



da 
w 



0. 



(A.31) 



When the Wilson loop average is saturated by the minimal surface [Eq. (|2.15j) ] the 
left-hand side of Eq. flA.30|) becomes proportional to the commutator of the functional 
derivatives applied on the minimal area. Using Eq. (|A.27[) . one obtains: 



6 



5 



5xP(a', 0) Sx a {a, 0) 5x a {a, 0) 5xP{a', 0) 



d_ 
da~ 



S(a' 



a 



da^(y)F(y-x(a,0)) 



_d 
da 



A(C) = 

-5{a'-a) J da^(y)F(y-x(a',0)) 



49 



+*(<7-<tV>,0)[ g j d<7°"(y)F(y - x(a, 0)) 

The right-hand side has the same structure as that of Eq. (jA.30|) and since the minimal 
surface satisfies the Bianchi identity (J3.10|) the same operations as above lead to the result 

(<fo/V,0) fe Q (a,0) ~ Sx a (a,0) 5x^(<r',0)) A ^ = °' (A ' 33) 

[Actually, it is the finite part of the quantity 9x 0® a ) / da au (y)F(y — x(a, 0)) that is relevant 
for the nontrivial content of the Bianchi identity, the singular part satisfying the Bianchi 
identity as an identity, independent of the type of the surface.] 

The diagonal constant metric (|A.4|) can also be used for the study of the forms of min- 
imal surfaces in simple cases. In that metric the equation of the minimal surface is given 
by the Laplace equation (|A.24[) . Choosing a plane as a reference surface, the orthogonal 
deformations of the surface with respect to the plane can be represented with functions of 
the type (jA.26|) . while the parallel deformations should be constructed according to the 
orthonormality conditions (jA.4|) . However, the sole knowledge of the orthogonal defor- 
mations is sufficient to have a rough idea of the form of the surface. We apply this study 
to one particular case of interest. 

Let us consider the case where the contour C of Fig ^ has been restricted in the 
following way: the line Xix' l x' 2 X2 represents now the three sides of a rectangle (lying in a 
plane) while the line x 2 X\, denoted C 2 i, is arbitrary but lying in an orthogonal plane to 
the rectangle and having in it the equation u(a, 0) = a(a) (Fig. |7J). 




x 2 




Figure 7: Contour C made of the three sides x\x' l x' 2 X2 of a rectangle and of the arbitrary 
line C21 lying in an orthogonal plane to the rectangle. 



Designating by u(a, r) the orthogonal deformation of the minimal surface with re- 
spect to the plane, its expression can be obtained from Eq. (jA.26|) by using the Fourier 
decomposition of the function a(a): a(cr) = Y^Li a n sin(?T,7ra). One obtains: 

00 sm(niia) sinh (mr(l — t)/A) 
u(a, t) = Y1 a « . r , V 7tt (A.34) 



sinh(mr/A) 
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In general this minimal surface does not have a simple form; u is different from zero all over 
the rectangle, but is mainly peaked in the region r ~ 0. Let us next take the limit x 2 — > oo, 
or A — > 0. Then u goes to zero for any r^O. Therefore, in that limit, the minimal surface 
shrinks to the surface made of the union of the rectangle xix' 1 x' 2 X2Xi and of the minimal 
surface bounded by the straight segment X2X1 and the line C21. Taking the time interval 
T proportional to Vi? in that limit, we see that T is now contained in the contribution of 
the rectangle only. Since the energy spectrum is provided by terms proportional to T, we 
conclude that the energy spectrum provided by the previous contour is the same as that 
of the rectangle. Therefore the shape of the contour C21 has no influence on the energy 
spectrum and contributes only to the wave functional of the corresponding bound state, 
made here of a static quark- ant iquark pair. The meaning of this result is that minimal 
surfaces, irrespective of the shapes of their contours, cannot produce quantum string- 
like excitations in the bound state energy spectrum, since rectangles do not contain such 
excitations. This was an expected result, for the minimal surface is the classical action of 
the open string and quantum excitations of it should only be searched for in fluctuations 
of surfaces about the classical trajectory |2H 1231 125] : the present derivation, however, 
provides another insight into the same property. 



B Normalization of the wave function 



In order to derive the normalization condition satisfied by the wave function, we start 
from Eq. (|P|) : 



Pl - (h w + h 



20, 



G 



1 5 n r 1 . 5 



+ 72L72 / daa 



6x0 {a] 



G 



-z<5 4 (xi - x'Jjil Goj - iS 4 (x 2 - x' 2 )i2L G ifi 



(B.l) 



i=l 



The series Y^i Gi y o and Yj Go,j that appear in the right-hand side are reminiscent of the 
series expansions of the quark and antiquark propagators in the presence of the gluon 
field (Sees. HJandEJ). They have the definitions: 



Y G ifi = -(tr c U(x2,x 1 )S 1 (x 1 ,x' 1 )U(x' 1 ,x 2 ))A, 
■■ (tT c U(x 2 , xi)U(x 1 ,x 2 )S- 2 (x' 2 , x 2 ))a- 



1=1 

00 



G o,j 
i=i 



:b.2) 



They can be interpreted as generalized gauge invariant quark (antiquark) propagators in 
the presence of the antiquark (the quark) at position X2 (xi)- The detailed study of the 
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structure of the series J2i Gifl an d J2j Go,j by means of the general formula ()5.6|) shows 
that they involve surface terms formed by the propagator lines and the external particle 
position, which implies that they contain, in addition to self-energy terms, interaction 
terms with the external particles through the functional derivatives of these surfaces. 
However, simplifications occur when one considers the equal time limits x 2 l = xil and 
X 2L = x il- Because of the presence of the four- dimensional delta-functions, the corre- 
sponding surfaces lie now in planes orthogonal to the timelike surface lying to infinity and 
containing X\X 2 (or are integrated around such planes). Such types of surface are generally 
negligible in front of surfaces that lie to infinity. In that case the main contributions that 
arise in the above series are those coming from the segments representing the boundaries 
of those surfaces and producing self-energy corrections. Those are calculated in Sec. 
[Eq. ()7.6j) ] and are found to be independent of the longitudinal momentum. We therefore 
approximate the right-hand side quantities of Eq. (jB.ljl by the quark propagators with 
self-energies. We have in momentum space 



J2G h0 ~-N c S( Pl , mi ) = -N c 



i=l 

00 



7i.pi - mi - S(pf,mi) : 



£ G 0tj ~ N c S(-p 2 , m 2 ) = N c — — -. (B.3) 

~{ -72-P2 -m 2 - £(-£>2 , ^2) 

The equal-time limit amounts to integrating with respect to the realtive energy variable 
in momentum space. Defining 

G ^ T ^ = lw)W) G(P ' p ' p,) {BA) 

(after having removed common factors (27r) 4 <5 4 (P — P')), and integrating Eq. ()B.1|) in 
momentum space with respect to pl and p' L we find: 

P L -h-h 2 - iK]G= l -N c ^ + ^ 7 iL72L(27r) 3 5 3 (/ - p' T ). (B.5) 

In the left-hand side we have designated by %K the kernel that results from the functional 
derivatives and which will survive in the bound state limit; it is essentially represented 
by the vector potentials calculated in Sec. El furthermore, we have extracted from the 
above functional derivative contributions, as in Sec. the self-energy parts and included 
them in the free Dirac hamiltonians hi and h 2 [Eqs. (j7.9j) - (j7.11j) ]. Next, we multiply both 
sides by G7il72l|(^- + |f) and take the bound state limit, represented in G by a pole 
in s = P 2 : _ 

G ~ -i-W-. (B.6) 

s^s S — 3q 
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Deriving both sides with respect to s and noticing that the left-hand side operator in Eq. 
(jB.5|) annihilates the bound state wave function [Eq. (J7.12J) ]. one obtains: 



d 3 p 7 



1 / hi h<2 



E-, 



^(P L ,p T ) i> 



2P L N c tp1pJ 1L ~f 2 L 



2 ^E\ E2 



7il72L- 



;b.7) 



(The trace is over the Dirac spinor indices.) To simplify the right-hand side, we notice 
that in the absence of the rotational motion of the flux tube, the vector potentials reduce 
to their longitudinal components [Eq. (J6.23J) ]. In that case, the wave function satisfies 
the equation 

Hpil) ~ e(p 2L )}^ = 0, (B.8) 
similar to the Salpeter equation case [SO], which implies that the right-hand side op- 

2 

7il72L, can be replaced by 1. Adopting this 



erator in Eq. (jB.7|) . 7il72l 
approximation, we obtain for the normalization condition: 



III 

E 2 



d 3 p 



T 



If h x 



(2^2^ 



E 2 



4> = 2P L N C 



;b.9) 



This formula is the same as that obtained from the Salpeter equation |3U| . 



C The Breit approximation 

The Breit approximation consists in transforming the wave equation (j7.12j) into a local 
equation in x-space. This is achieved by first replacing the energy sign operators appearing 
in the potentials (Jfi.20)) by 1; in the nonrelativistic limit this approximation affects only 
terms of order 0(l/c 4 ) or higher. Second, the operators E a , a = 1,2, and L 2 /x 2 that 
appear in the potentials are replaced by their mean values. Third, the self-energy terms 
are neglected. The wave equation takes then the form ()6.13|) . where h a0 , a = 1,2, are the 
free Dirac hamiltonians [Eqs. ()6.7j) ] and potentials A afl , a = 1,2, have been submitted to 
the approximations described above. In order not to spoil the property of spontaneous 
chiral symmetry breakdown realized by the self-energy parts of the Dirac hamiltonians 
h a (a = 1,2) [Eqs. ()7.9|) ] in the massless limit, we replace the total cm. energy P L 
by \J P 2 + 4:Eq, where E is a constant the value of which is fixed by the requirement 
that the ground state be massless. Those approximations do not affect the high energy 
properties of the spectrum, since these are mostly determined by semi-classical properties 
(in particular the Regge behavior). 

The resulting equation is now a Breit type local equation in x that can be solved with 
standard methods. In this form, however, a difficulty arises coming from the fact that 
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vectorlike local potentials of the confining type at the classical level do not seem to confine 
at the quantum level due to instabilities that appear at large distances [23] • The difficulty, 
however, is only apparent. The study of the Pauli-Schrodinger equation that results from 
the Breit type wave equation shows that the Darwin terms of the equation contain a 
repulsive singularity at a finite distance, depending on the energy of the bound state, which 
confines the quarks inside a baglike system with an energy dependent radius and forbids 
any tunnelling between the inner domain of the bag and its outer domain. Therefore, 
confining normalizable solutions do exist; they are obtained with the boundary condition 
that the component ip3 of the wave function [Eq. (J8.14)) ] vanishes at the boundary of the 
bag. That boundary condition, breaks, however, implicitly chiral symmetry, since it does 
not operate in a symmetric way between chiral partners. This is why, in spite of the facts 
that the whole interaction is vectorlike and self-energy parts have been neglected, the 
spectrum displays chiral asymmetry and a massless ground state appears if the constant 
E is correctly adjusted. 

The existence of that type of confining solutions with vectorlike potentials was already 
pointed out by Geffen and Suura jHS]. In that work, the masslessness of the ground state is 
ensured by the presence of the short-distance Coulomb-like potential and the adjustment 
of the corresponding coupling constant. 
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